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^ ^ We present a self-contained summary of the theory of hnear cosmological per- 

^ ^ turbations. We emphasize the effect of the six parameters of the minimal cos- 

I mological model, first, on the spectrum of Cosmic Microwave Background tem- 

Q perature anisotropies, and second, on the linear matter power spectrum. We 

^ briefly review at the end the possible impact of a few non-minimal dark matter 

^ and dark energy models. 
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Purposes 

Cosmology is progressing by leaps and bounds thanks to a spectacular 
amount of observational data. It provides crucial clues for particle physics, 

fS| and more generally for high energy physics. A key role is played by Cosmic 

Microwave Background (CMB) and Large Scale Structure (LSS) data. Their 
interpretation relies on the theory of cosmological perturbations. 
• • The aim of this course is to introduce this theory, staying at the linear 

. ^ level, and hopefully with a simple, self-contained and original approach. It 

is targeted at PhD students in astrophysics and cosmology, as well as re- 
^ searchers in particle physics willing to follow current developments in cos- 

mology. We discuss in particular the effect of the parameters of the minimal 
cosmological model on the spectrum CMB temperature anisotropies, and 
on the matter power spectrum. We briefly review at the end the possible 
impact of non- minimal dark matter and dark energy models. 

These notes are based on a cycle of five lectures of 90 minutes. Since 
this is not much for introducing linear cosmological perturbation theory, 
the lectures remained at a qualitative level, with more emphasis on physical 
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ideas than on the algebra. These written notes contain roughly the same 
amount of information as the lectures. After reading these few pages, the 
reader may wish to improve his/her level of understanding, while sticking 
to the same structure and same notations. For this purpose, we suggest to 
read chapters 2, 5 and 6 of Ref. 1: they contain a similar presentation, still 
simple and qualitative, but more detailed. Optionally, the reader can also 
grab there a description of neutrino effects, which are completely neglected 
here. 

For a quantitative description of perturbation theory, gauge transforma- 
tions, CMB physics and linear or non-linear structure formation, the reader 
should refer to the specialized literature. In the bibliography, we list a very 
incomplete selection of outstanding references. 

Sections 2 and 3 of these notes are illustrated with several figures that 
represent the qualitative evolution of cosmological perturbations. It would 
have been easy to produce them with a Boltzmann code, as we did in Ref. 1. 
We deliberately chose instead to include only scanned hand- drawings. The 
first reason is that drawings offer the possibility to exaggerate some effects 
for pedagogical purposes: they are often more illuminating than exact nu- 
merical plots. The second reason is that they better render the conviviality 
of the blackboard presentation at TASI. 

1. Linear cosmological perturbations 
1.1. Classification 

We decompose the metric and stress-energy tensor of the universe into 
spatial averages and linear perturbations. 



where g^jy stands for the metric of the homogeneous and isotropic 
Friedmann-Lemaitre (FL) model. Being symmetric, the two perturbed ten- 
sors contain ten degrees of freedom each, describing different aspects of 
gravity. Bardeen showed in 1980 that they can be decomposed on the ba- 
sis of scalars, vectors and tensors under spatial rotations (spatial rotations 
play a special role because they leave the FL background invariant). These 
three sectors are decoupled at first order in perturbation theory. 

In the vacuum, scalar and vector perturbations vanish, while tensor per- 
turbations can propagate if they have been excited: they account for gravita- 
tional waves, the only "real" (propagating) gravitational degrees of freedom 



(1) 
(2) 
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in General Relativity (GR). In the presence of matter, scalars represent the 
response of the metric to an irrotational distribution of matter, and gener- 
alize Newton's theory of gravitation. Vectors represent the response of the 
metric to vorticity, and describe phenomena with no equivalent in Newton's 
theory, called "gravito-magnetism" . 

In minimal cosmological models, the vorticity of the various matter com- 
ponents decays with time, and vectors can be neglected. Tensors may play 
a small role in CMB anisotropics, that we will mention briefly in Sec. 2.8. 
They can be studied separately, since they decouple from the scalar sector 
at first order in perturbations. Hence this course will be essentially focused 
on scalar perturbations. 

The four scalar components of both the metric and stress-energy per- 
turbed tensors are contained in: 

(1) the (00) term, 

(2) the trace of the {ij) matrix, 

(3) the irrotational part of the (Oi) vector, 

(4) the traceless longitudinal part of the (ij) tensor. 

For the perturbed metric Sg^i,^ these components correspond (in the same 
order) to: 

(1) the generalized gravitational potential t/^, 

(2) the local distortion (j) of the average scale factor a{t): the "local scale 
factor" is given by (1 — (/))a, 

(3) the potential b such that Sgoi = dib, 

(4) the potential /i of the metric shear: Sgij = {didj — \5ij/S)ii . 

For the perturbed stress-energy tensor dT^j^^ these components (still in the 
same order) represent: 

(1) the energy density perturbations dp (we will usually refer to the relative 
perturbations 5 = 5p/p)^ 

(2) the pressure perturbations 

(3) the potential v of the irrotational component of the flux of energy, 
6T^ = diV {v is sometimes called the velocity potential, since in the 
case of a fluid it is related to the bulk velocity), 

(4) the potential s of the shear stress or anisotropic stress: ^TJ = {didj — 

It is equivalent to use as a variable the "velocity potential" v or the "velocity 
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divergence" 6 defined as 

di5if = {p+p)e = ^v. (3) 

Similarly, we can use the function a instead of s, with the definition 

{didj - A)(5T; = {p+ p)a = A(As) . (4) 

The function a is usually called the anisotropic stress, although the true 
anisotropic stress is the component of STj deriving from the potential s. 
The factors {p-\-p) in the two previous equations are introduced in the def- 
initions in order to obtain simple equations (some authors use alternative 
notations without these factors or with different ones). To summarize, we 
see that we can manipulate four degrees of freedom representing the scalar 
perturbations of matter fields, that can be chosen to be the density fluc- 
tuation, pressure perturbation, velocity divergence and anisotropic stress: 
{S,Sp,0,a}. 

1.2. Gauges 

In an idealised FL universe, there is only one time slicing compatible with 
the assumption of homogeneity. Instead, in a perturbed universe, there is 
an infinity of time slicings compatible with perturbation theory (i.e. such 
that on each slice, all quantities remain close to their average value). 

The perturbation of any quantity in a given point is the difference be- 
tween the true and the average quantity in this point. For instance, for the 
total energy density p, 

Sp{t,x)=p{t,x)-p{t) . (5) 

While p(t, x) is a locally, unambiguously defined quantity, p{t) depends on 
the choice of equal-time hypersurface going through the point (t, x). With 
a different choice, p{t, x) would be compared to the average performed on a 
different sheet, that would take a different value. Hence Sp{t, x) also depends 
on the choice of time slicing. 

A gauge is a choice of time slicing. Gauge transformations are induced 
by coordinate transformations x^-\- mapping the points of one time 

slicing to those of another time slicing. All coordinate transformations do 
not induce a valid gauge transformation: the condition that perturbations 
must still be linear after the transformation restricts to be small in every 
point. 

A naive study of the equations of motion of perturbed quantities would 
be plagued by the freedom to change the gauge without changing physical 
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results: some solutions of the full equations would be "gauge modes" with 
no observable consequences. To deal with this issue, one can adopt one of 
two point of views: 

• one can derive gauge- invariant quantities (i.e., non-trivial integro- 
differential combinations of the metric and stress-energy tensor com- 
ponents left invariant by a gauge transformation), and gauge-invariant 
equations of motions for these quantities. Note that there are four scalar 
degrees of freedom in ^^^^^ and two scalar degrees of freedom in the 
four- vector field inducing gauge transformations: namely, eo and the 
potential e such that ci = d^e. Hence we can use gauge transformations 
to cancel two scalar degrees of freedom, and build up two indepen- 
dent gauge-invariant scalars. One way to define them is through the 
two Bardeen potentials $a and defined by Bardeen (1980) as two 
integro-differential combinations of (j), fi. 

• one can fix the gauge, i.e. introduce a condition such that the time 
slicing is unique. Then the number of independent solutions to the 
equations will be the same as in the gauge- invariant formalism. In any 
case, one can show that truly observable quantities are always inde- 
pendent of the gauge. Obtaining them after solving equations in the 
gauge-invariant formalism or in one particular gauge should not make 
any diff'erence in practice. 

A convenient gauge choice for a pedagogical introduction to CMB and 
matter scalar perturbations is the so-called Newtonian gauge or longitudinal 
gauge, in which one imposes that non-diagonal scalar perturbations of the 
metric vanish: b = fi = 0. This prescription can be showed to fix a unique 
time slicing. In this gauge, adopting units such that c = 1 and using proper 
time t, the line element reads: 

ds'^ = -(1 + 2i;)dt'^ -h (1 - 2(l))a^dP , (6) 

where dl ^ stands for the cartesian measure (dx^ + dip' + dz^) for a fiat FL 
model, or for (1 zb kr^Y^dr^ + r^(sin^^d^^ + d^p^) for an open/closed FL 
model in spherical coordinates. We are still free to redefine time (by defini- 
tion, a time redefinition leaves the time slicing invariant). Lots of results in 
cosmological perturbation theory look simpler when using conformal time 
7/, defined up to a constant by dr\ = dtja. In this course we fix the con- 
stant in such way that 77 — > at the vicinity of the initial singularity, when 
p — > 00 (this prescription would not work if we were studying cosmological 
infiation, but in this course, we are not). Conformal time is convenient be- 
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cause photons traveling in a fiat unperturbed FL universe along geodesies 
crossing the origin of the system of coordinates obey to dr = dr] (this 
comes from adl = dt, i.e. from ds = with the restriction dO = dc/) = 0). 
Hence conformal time is a measure of time based on the comoving dis- 
tance travelled by a photon^, and the comoving distance to a given object 
coincides with its "look-back conformal time". In this course we will use 
dots for derivatives with respect to proper time and primes for derivatives 
with respect to conformal time. The Hubble parameter (or expansion rate 
parameter) reads 

H='- = ^,, (7) 

the Hubble radius is Rh = and the condition that a Fourier mode of 

physical wavelength A crosses the Hubble radius is 

27r 1 , a' 

A = Rh ^ = — ^k^aH=— . (8) 

k H a 

One advantage of the Newtonian gauge is that the gauge- invariant 
Bardeen potentials {^a^^h} reduce in this gauge to the metric pertur- 
bations {(j),2p}: the evolution of the Newtonian metric perturbations in- 
forms us directly on that of two gauge- invariant quantities. Other inter- 
esting properties of ip} appear when writing the Einstein equations in 
the Newtonian gauge. The full Einstein equations linearized at first order 
in perturbations feature four equations relating scalar degrees of freedom. 
One of them, associated to the traceless longitudinal part of SGj = SttGTJ, 
gives (assuming a flat FL background) 

?(fc/a)2(^ -^P)= STrG^fe +p,)a, , (9) 

X 

where the index x runs over all the species contributing to the total stress- 
energy tensor. This means that when the universe contains only shearless 
components with cix = (as would be the case in the presence of perfect 
fluids), the two metric perturbations are equal. Next, the Einstein equation 
SGq = SttGTq gives (still assuming a fiat FL background) 



2a" 



a 



8^GVp,(5, . (10) 



The term on the right-hand side involves the total energy perturbation 
^Ptot = J2x Px^x- In the short scale (more precisely, sub-Hubble) limit, the 



^In a universe with non-zero spatial curvature, this remains true, provided that the 
comoving distance is defined not hke r, but hke x = /(I ^ kr'^)~^^'^dr. 
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term containing k'^ dominates the other terms in the square brackets, and 
we recover the Poisson equation 

-^.^ = 47rG5ptot , (11) 

where the factor —k'^/a^ represents the Fourier transform of the physical 
Laplace operator in an expanding universe. Note that in the Poisson equa- 
tion one may have expected to see the generalized gravitational potential 
^p instead of 0: however, in the sub-Hubble limit, the shear of individual 
components is usually either null or negligible, so that (j) = xp. 

1.3. Equations of motion 

In the minimal cosmological scenario, the universe features several species 
with spatial fluctuations, described with different equations because of their 
distinct properties: cold dark matter (CDM) is non-relativistic and collision- 
less, neutrinos are ultra-relativistic and collisionless at the times of inter- 
est, baryons are non-relativistic and smoothly interpolating from a strongly 
coupled to decoupled regime, and finally photons are ultra-relativistic and 
interpolating between the same two regimes. 

The equation of conservation of the total stress-energy tensor, DfjTj^ = 
(deriving from Bianchi identities), yields two scalar and two vector equa- 
tions. The scalar ones are the conservation of energy equation and the Euler 
equation. 

For an single component experiencing no interaction with other species 
(other than gravitational), the equation D^Tj^ = applies to the individual 
stress-energy tensor: it gives one continuity and one Euler equation for that 
component. The evolution of single component experiencing interactions is 
also given by the continuity and Euler equation, but with an extra source 
term accounting for stress-energy injection/leak caused by the interaction. 

We have seen that the perturbations of each component x can be de- 
scribed by four variables {d^^ Spx, ^x^^Fx}- Hence, in general, two equations 
of motion are not sufficient for closing the system. However: 

• for a perfect fluid, microscopic interactions impose local thermodynam- 
ical equilibrium. The pressure is then isotropic^, with ax = 0. In ad- 



^An anisotropic shear stress a ^ reflects the fact that in each given point, particles 
travel with different velocities (due to some intrinsic velocity dispersion and/or a super- 
position of several flows in phase space), leading to anisotropic pressure. This contradicts 
the assumption of a perfect fluid, in which local interactions result in a unique bulk ve- 
locity (after coarse-graining over microscopic scales), and erase anisotropic pressure. 
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dition, pressure perturbations obey to 5px = c^dpx^ where Ca is the 
adiabatic sound speed inferred from the equation of state of the fluid. 
If ax vanishes and 6px is a function of Spx^ perturbations in the fluid 
are described by only two independent functions 6x = Spx/px and Ox- 
If the collision term also vanishes or is specified, the two equations of 
motion (continuity and Euler) are sufficient for closing the system and 
computing the evolution of perturbations. 
• for a decoupled or weakly interacting species, there are no such simpli- 
fications concerning the anisotropic stress and pressure perturbation. 
Hence, in general, the two equations inferred from stress-energy con- 
servation are not sufficient. For such species, one has to use the more 
general Boltzmann equation, giving the evolution of each phase-space 
distribution: 



where the sum holds over the species y coupled with x. Each phase- 
space distribution can be decomposed into a background and pertur- 
bation part: 



where p stands for momentum, \p\ for its modulus, and the background 
part does not depend on the direction of p by assumption of isotropy. 
• fully decoupled CDM is a particular case of a collisionless species with 
negligible velocity dispersion (the word "cold" refers precisely to this 
last assumption). Hence, it behaves in the same way as a pressure- 
less perfect fluid, although in reality it has no interactions and should 
not be called a fluid. Since the velocity dispersion is negligible, in a 
given point, all particles share the same velocity, imposed by gravita- 
tional flows (while for non-cold collisionless species, the velocity would 
get two contributions, one from gravitational flows, and one from the 
phase-space distribution function) . Hence the anisotropic pressure van- 
ishes (see the previous footnote) . The pressure perturbation Spx is also 
related to the velocity dispersion in a given point, anf can be neglected 
with respect to the density perturbation in the CDM case. Hence CDM 
is formally equivalent to a perfect fluid with no anisotropic stress and 
no pressure perturbation (or in other words, with a sound speed Cg = 0). 
In that case, the two equations of motion inferred from D^jT}^ = are 
sufficient, like for a fluid. 




(12) 



fx iv, x,?) = fxivApl) + ^fx iv, ^, p) 



(13) 
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Gravitational interactions between species are accounted by the presence of 
metric perturbations in each equation of motion (more specificahy, terms 
in k'^ip accounting for gravitational forces, and terms in (p' accounting for 
dilation effects, i.e. for local distortions of the scale factor with respect 
to a). Hence, in order to close the full system of equations, we still need 
two independent relations, to be chosen among the four scalar Einstein 
equations: they provide the value of k'^ip and 0' at each time, as a function 
of all matter fields. 

1.4. Initial conditions 

We wish to study the evolution of matter perturbations, starting from some 
early time at which all Fourier modes of interest (those which are observ- 
able in the CMB spectrum and in the matter power spectrum on linear or 
mildly non-scales) are still outside the Hubble radius. Indeed, super-Hubble 
modes experience a trivial evolution, unaffected by small scale interactions 
(Thomson or Coulomb scattering, usual gravitational force V0, etc.) Hence, 
the perturbations evaluated at some arbitrary time but on super-Hubble 
scales reflect directly the mechanism responsible for the formation of per- 
turbations in the very early universe. In the standard cosmological model, 
these initial conditions can be inferred from inflation. 

Typically, a good time for setting initial conditions is when the redshift 
z = ao/a — 1 is of the order of 10^: at this time, all comoving scales that 
are observable in the CMB and linear matter power spectrum still verify 
k < aH. 

It is crucial to understand that, as long as the background cosmology 
is assumed to be of the FL type, the perturbed stress-energy momentum 
tensor ^T^ must be diagonal on super-Hubble scales. Indeed, the back- 
ground tensor is diagonal, and of the form: diag(— p,^,^,^). This can 
be showed to be the most general assumption compatible with homogeneity 
and isotropy. Let us assume that we Taylor-expand (5T^ in powers of the 
variable {kr]). For any power-law scale factor, aH is given by 1/?^ times a 
factor of order one. Hence the limit kr] <^ 1 represents precisely the super- 
Hubble limit. In the Taylor expansion, the zero mode should share the same 
properties as the background solution, and be diagonal. Higher order terms 
account for contributions to STj^ growing with time, and possibly becoming 
important around the time of Hubble crossing. 

The total scalar perturbations Sp and Sp are the only one preserving the 
diagonal form of STj^: we conclude that they are the only one that do not 
vanish at order zero in the {kr]) expansion. Using stress-energy conservation 
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equations, one can show that the part of STj^ associated with the velocity 
divergence is of order one in (/cry), while the part associated to the total 
anisotropic stress is of order two. 

Suppose that the universe contains initially N uncoupled perfect fluids^, 
with known sound speeds = dp^/^Px- There are 2N independent 
initial conditions, corresponding to possible initial values of each 5^ and 
each 5'^. Importantly, in the 2 TV- dimensional basis of IC's, one basis vector 
is very special, as we shall see below. 

Before studying perturbations, one should have solved for the back- 
ground evolution: all background quantities should be known, including for 
instance the density px{vi) and pressure pxij]) of ^ach species x. Now, let 
us assume that the real universe is perturbed initially by a single degree of 
freedom (one may say, by a single initial time shifting function) . This is the 
case in single-field inflationary cosmology: during inflation, there is a sin- 
gle clock (the infiaton), and perturbations arise from a single time shifting 
function (the infiaton perturbation). 

As long as we are dealing with super-Hubble modes, we can neglect 
microscopic interactions and say that the evolution in each point (in fact, 
in each Hubble patch) is still given by homogeneous cosmology, taking this 
shift function 5iri{x) into account: 

Vx, px{r], x) = p^{r]^ 5r]{x)) = p^iv) + Pxiv) ^v{^) 

Px{r],x) =Px{r] + 5r]{x)) = Px{v) ^ Pxiv) ^vi^) (14) 

where in the last equalities, terms of order two or higher in 6r] have been 
neglected. The above ansatz restricts a lot the choice of possible initial 
conditions. Indeed, for each uncoupled species, the energy conservation law 
for the background gives 

p', = -3"^{p,+p,) . (15) 
Combining this equation with the above ansatz, we get 

v., = -3^% = -3^ Srjix) . (16) 

Px-\-Px Ci Px 

Since the last term is independent of the index x, we see that for each pair 



^The discussion presented in this section could be generalized to N coupled species, not 
all of them being perfect fluids: the conclusions would not change qualitatively, and we 
restrict here to N uncoupled fluid for simplicity. 
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of fluids, 

Px -^Px Py ^ Py 

Note that this famous relation is a consequence of Eqs. (14). It shows that in 
presence of such initial conditions, everything is flxed up to a single function 
of X. Let us take the example of a universe containing only photons, baryons, 
cold dark matter and neutrinos. We can use the fact that for non-relativistic 
species Px ^ Px^ while for ultra-relativistic ones px = Px/^- Hence, if one 
function is known — for instance, Sj{x) at initial time — the others can be 
derived from 

3 3 

Sh = ^cdm = -^Sr, = -S^ . (18) 

It is a simple exercise to prove that the ansatz of Eqs. (14) implies additional 
important relations: 

• flrst, 

yx, 5pM x) = Sp:,{v, f) , (19) 

Px \V ) 

i.e., as long as we are dealing with super-Hubble modes and Eqs. (14) 
are satisfled, each component features an adiabatic sound speed inde- 
pendent of X, obeying to 4(r/) = p^M/pxiv)^ 

• second, 

^Ptot {r],x) = cl {f]) 6ptot {r],x) , (20) 

where the squared total sound speed Cg can be easily expressed as a 
weighted average over the squared adiabatic sound speed of individual 
components. Hence, the total matter content resulting from the sum of 
all components also features an adiabatic sound speed. 

For the reason expressed in the last item, such initial conditions are usually 
called adiabatic initial conditions. The word adiabatic refers to thermody- 
namics, but it is actually more helpful to think of adiabatic initial conditions 
as resulting from the assumption of a single inhomogeneous degree of free- 
dom in the early universe: then Eqs. (14) stands for the true deflnition of 
adiabatic initial conditions, and Eqs. (17, 18, 19, 20) for its consequences. 

More general initial conditions not obeying to Eqs. (14) can be expanded 
on diff"erent bases. A famous basis is formed by (i) the adiabatic mode; (ii) 
{N — 1) non-decaying isocurvature modes, getting their name from the 
property that for each of them, the total density (and spatial curvature) 
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perturbations vanish asymptotically in the super-Hubble limit, while two 
species have opposite density perturbations compensating each other; (iii) 
TV decaying modes that are irrelevant for most purposes. 

It is clear from the previous discussion that non-adiabatic initial con- 
ditions should only be considered when assuming that primordial pertur- 
bations are generated by more than one degree of freedom (for instance, 
two inflaton fields, or one infiaton and one axion, etc.) The assumption 
of several degrees of freedom is necessary but not sufficient. The primor- 
dial universe may contain a mixture of adiabatic and isocurvature modes 
until a time at which all species are brought in thermal equilibrium. At 
that time, if we further assume that all chemical potentials vanish, the 
perturbations of each species can be inferred from those of temperature, 
T(ri^x) = T{rf) -h 6T{r],x). Then, temperature plays precisely the role of a 
single time-shifting function. Any non-adiabatic initial condition is washed 
out and becomes irrelevant. Hence, isocurvature modes can be observable 
— for instance in the CMB — only under additional assumptions. For in- 
stance, one species carrying isocurvature perturbations might remain out of 
equilibrium at all times, or might feature a chemical potential with spatial 
fluctuations. There exist a few non-minimal, but still reasonable scenarii fea- 
turing isocurvature perturbations (with axions, curvatons etc.). Ultimately, 
the presence of isocurvature modes is to be tested with observations. Cur- 
rent CMB observations put strong limits on the amplitude of such modes, 
and prefer purely adiabatic initial conditions. We will restrict to this case 
in the rest of this course. 

For adiabatic initial conditions, we found the relation (17) holding be- 
tween density fluctuations. Also, if we do not consider the case of species 
with a non-zero anisotropic stress, we can assume that = -0 at initial 
time. By plugging Eqs. (18), (20) and (/) = ip into the Einstein equations, 
one is led to a second-order differential equation for tjj only. During the ra- 
diation dominated era, one can show that this equation has two solutions, 
one constant in time, and one decaying. The (00) Einstein equation then 
shows that the constant solution is related to density fluctuations through 

—2ip = —2(j) = 6tot 6j = constant . (21) 

Hence, in the Newtonian gauge and for adiabatic initial conditions, super- 
Hubble metric fluctuations and density fluctuations are static. In fact, on 
super-Hubble scales, there can only be some time evolution when the uni- 
verse changes of total equation of state (or equivalently, of expansion law). 
This is the case at the time of equality between radiation and matter. Dur- 
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ing matter domination and on super-Hubble scales, they freeze out again. 
Then, the relation — 2-0 = —20 = ^tot and Eq. (18) are still satisfied, but 
now Stot — = ^cdm = |^7- This detail will become important when 
studying the Sachs- Wolfe effect in section 2.3. 

1.5. Power spectra and transfer functions 

The theory of cosmological perturbations is a stochastic theory: the fluctu- 
ations of a given quantity in a given point, A{r],x), obey to a distribution 
of probability. As long as we stick to linear perturbation theory, there is 
a "linear transport of probability" from one time to another time. Let the 
probability of ^4 in a given point at time r]i be Vi{A). In the same point 
and at 772, the linear evolution would have transformed each value A into 
aA, where a is the linear growth (or decrease) factor of A between r]i and 
7/2- So the probability of A at time rj2 is given by V2{A) = Vi{A/a). This 
"linear transport of probability" implies that the linear evolution respects 
the shape of the probability distribution, and rescales all its statistical mo- 
ments of order n by a^. In particular, if the initial probability is Gaussian, 
the statistics will remain Gaussian at any later time, and the evolution of 
the system can be formulated as an evolution of the root mean square of 
all quantities. In summary, as long as we assume linear perturbations with 
Gaussian initial conditions, our goal is to solve for the evolution of the root 
mean squares of the fluctuations. 

The equal-time 2-point correlation function of any quantity A in real 
space is given by 

{A{rj, x)A{rj, f )) = ^V, x, f ) ^V, W - ^1) , (22) 

where the last equality holds as a consequence of Statistical Homogeneity 
and Isotropy (SHI, assumed to hold in a perturbed FL universe). Indeed, 
the correlation function should be invariant under spatial translations and 
rotations. 

We can go to Fourier space and use the same letter to denote the Fourier 
transform of A. If A is real, A{r], —k) = ^*(77, k). It is easy to show that the 
previous equation (using the assumption of SHI) implies that the two-point 
correlation function in Fourier space reads 

{A{r], k)A*{r], k')) 5n{k' - k)PA{k) , (23) 

where Sd is the Dirac distribution. Here, statistical homogeneity is respon- 
sible for the fact that the two-point correlation vanishes for k ^ k' ^ and sta- 
tistical isotropy for the fact that Pa only depends on the modulus k = |A:|. 
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The function Pa(^) is usually called the power spectrum of A. Some authors 
prefer to refer to the "dimensionless power spectrum" , defined as 

PA(fc) = ^PA{k) . (24) 

The reason is that typical expressions for the average of various quantitites 
in real space consist in the convolution of the power spectrum with some 
window function /(/c), like in 

j-^^p^(k)m = J -^p^(k)m = J diogk rA{k)m . (25) 

Hence the dimensionless spectrum VA{k) stands for the weight of each 
logarithmic interval in the integral. The term "scale- invariant spectrum" 
refers to 7^a(^) being independent of i.e. Pa(^) oc . 

We know that for adiabatic initial conditions, all perturbations are 
related to each other through Eqs. (18, 21). Hence, with Gaussian adia- 
batic initial conditions, specifying the primordial power spectrum of one 
quantity is sufficient for knowing everything about the system. For in- 
stance, if we assume that the power spectrum of the metric perturbation 

is a given function P^(/c), then we infer from Eqs. (18, 21) that the 
photon and baryon primordial spectra are given by Pj{k) = 4P^(/c) and 

Pb(fc) = ^p^{k) = ip^{k). 

By convention, the primordial spectrum is usually given for the variable 
7^, which represents the spatial curvature perturbation on one initial co- 
moving hypersurface (i.e. an hypersurface orthogonal to the energy flux of 
the total cosmic fluid in each point). The advantage of using this quantity 
is that it is conserved on super-Hubble scales for adiabatic initial condi- 
tions (while (j) and ijj get rescaled when the equation of state of the universe 
changes, e.g. at radiation-matter equality). In the Newtonian gauge, IZ re- 
lates to ijj and to the total density perturbation through 

n = ^-\^^^. (26) 

3 ptot + Ptot 

The power spectrum of a given quantity at some arbitrary time can be 
decomposed into two parts, one accounting for initial conditions, and one 
accounting for linear evolution with time: 



{A{7J,k)A*{7J,k'))=5D{k' -k) 



■R{k) 



Pn(k) . (27) 



In the above equation, there is no time argument for 7l{k) since this quan- 
tity is conserved on super-Hubble scales. We just assume that IZ is evaluated 
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at a time such that k <C aH for ah modes of interest. In a FL universe, 
the equations of motion of ah perturbations respect isotropy, and do not 
depend on the direction of the wavector k. Hence the ratio between brack- 
ets in the last equation is a function of k, not k. This function accounts 
for the hnear evolution, independently of initial conditions. It is called the 
"transfer function" of A. In this course, we will denote transfer functions 
with the same letter as the perturbations themselves, but with the modulus 
of k as an argument, i.e. 

In summary of this section, solving for cosmological perturbations (in a 
model with adiabatic Gaussian IC's) amounts in 

• postulating a primordial spectrum, or calculating it within the frame- 
work of a model, for instance of inflation; 

• solving the equations of motion of all perturbations, with quantities 
normalized initially to lZ{k) = 1, in such a way to obtain all transfer 
functions A{r], k) and to infer the evolution of the various root mean 
squares. 

2. CMB temperature anisotropies 

From now on, we will assume for simplicity that the universe is flat. 



2.1. Photon scattering rate 

CMB physics consists is the study of electron, baryon and photon per- 
turbations on cosmological scales, taking into account their gravitational 
coupling with collisionless species such as decoupled neutrinos and CDM. 
Electrons and baryons carry opposite electric charges and are coupled to 
each other through very eflicient Coulomb scattering processes. Electrons 
and photons are coupled through Thomson scattering, which is the limit of 
Compton scattering when electrons are non-relativistic and photons carry 
a smaller energy than the rest mass of the electron. Then, the scattering 
process results mainly in a deflection of the photon, with a negligible trans- 
fer of energy between the two particles. The leading interaction between 
baryons and photons is the gravitational one. 

In units such that c = 1, the Thomson scattering rate (with respect 
to conformal time) is given by F = cFTdrieXe^ where cft is the Thomson 
scattering rate, a is the scale factor, Ue is the total electron number density 
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(scaling like due to dilution), and Xe is the ionized electron fraction. 
The product aUe scales like a~^. The ionized fraction is close to one at high 
energy. Then, at the time of recombination between electrons and nuclei 
(around z ~ 1080), which falls at the beginning of matter domination, Ue 
drops abruptly to very small values. This causes Thomson scattering to 
become suddenly very inefficient, and photons to decouple from electrons. 

Hence photon decoupling is the story of Thomson scattering becoming 
inefficient, while Coulomb scattering remains very strong. For that reason, 
one can describe electrons and baryons as a single tightly-coupled fluid. 
People often refer only to baryons for simplicity. At early time, the full 
system of (electrons) -baryons-photons is also tightly coupled, but later on, 
it splits progressively into two collisonless species, (electrons)-baryons on 
the one hand, photons on the other hand. 

The thermo dynamical description of recombination is very technical, 
due to the different energy level of atoms (in particular, of hydrogen). In 
order to understand CMB anisotropics, we only need to describe the main 
results of recombination studies at a very qualitative level. 

• The free electron fraction Xe{r]) starts from one at high redshift. It 
decreases sharply at the recombination time (corresponding to z ^ 1080 
or T 0.3 eV), and freezes at a very small value (due to departure 
from thermal equilibrium). Around z ^ 10, star formation causes a 
reionization of the universe, and Xe goes up again to one. 

• The Thomson scattering rate F = aTCtUeXe evolves like a~^Xe. Before 
recombination, V ^ and the universe is opaque. The sudden drop 
of Xe at recombination renders the universe transparent: F <C ^. Due 
to the dilution factor coming from rig, F remains much smaller than ^ 
even during reionization, and the universe keeps being transparent (this 
is why despite of reionization, most photons emitted at recombination 
do not interact anymore, and allow us to observe anisotropics on the 
last scattering surface). 

• The optical depth r(7y) = diniV{rf) represents the opacity of the uni- 
verse at a given time, when seen from today (when 77 = r/o). It tends 
to infinity when rj — > 0, falls below one at recombination and stabilize 
at a value of the order of 0.1 between recombination and reionization. 
After reionization it decreases smoothly and reaches zero today by def- 
inition. 

• The visibility function g{r]) = —T'e~^ gives the probability that a CMB 
photon seen today experienced its last scattering at time 77. It starts 
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from negligible values at high redshift (suppressed by the factor). It 
has a narrow spike around the time of recombination, and then it falls 
again to negligible values due to the smallness of r' between recombi- 
nation and reionization. It develops a second smaller and wider spike 
around reionization. This function shows that most CMB photons did 
not interact between the last scattering surface and today, while a mi- 
nority rescattered at reionization. The width of the recombination spike 
gives an indication on the thickness of the last scattering "surface" . 
The diffusion length Ad(^) is an important quantity for understanding 
the damping of temperature anisotropics on small scales. At any given 
time, the comoving mean free-path (mfp) of photons is given by rmfp = 
l/Vij]) (still in units where c = 1). If the photons experience a random 
walk analogous to Brownian motion in a gas, the comoving distance 
over which they travel between time Tyini and 77 can be approximated 
by 

1/2 



J riini 



df] r r 



mfp 



nl/2 






/ df] T-^ 




-J mm 



(29) 



This physical diffusion length of photons is given in this approximation 
by Ad — ar^. It grows quickly from very small to very large scales 
(comparable to the Hubble scale) around the time of recombination. 



2.2. Boltzmann equation 

Since photons decouple from baryons near the recombination time, we can- 
not describe them with fluid equations, and need to solve the Boltzmann 
equation 

|^/ = C[/,/e] (30) 

at order one in perturbations. The right-hand side stands for the photon- 
electron coupling due to Thomson scattering. As explained in the previous 
section, electrons and baryons are so tightly coupled that it makes no differ- 
ence to think of this term as a photon-electron or photon-baryon coupling 
term. Solving this equation is involved because the photon phase-space dis- 
tribution / involves many arguments, /(r^, x,p). Fortunately one can reduce 
the dimensionality of the problem. First, we notice that as long as photons 
are in thermal equilibrium with electrons (and hence with baryons), they 
are entirely described in any point by the local value of the equilibrium tem- 
perature T(7y, x). The phase space distribution is then of the Bose-Einstein 
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form: 

f{r],x,p)= . (31) 

It can be expanded into a background part and a first-order perturbation, 
f = f^Sf, with: 

= , Sfiv^x^p) = - A . (32) 

e^M - 1 dlogp T[r]) 

We see that in the tightly-coupled regime we could replace the variable 
f{r],x,p) by the lower- dimensional variable ©(77, x) = [ST{r],x)/T{r])]. The 
Boltzmann equation leads to an equation of motion for 6(77, x). 

At later times, when photons decouple, the shape of the Bose-Eisntein 
distribution of photons is preserved. This can be infered from the geodesic 
equation. This equation tells how the individual momentum p of photons 
evolve when they are decoupled and they travel in the perturbed universe. 
In the Newtonian gauge, it reads 

^M = _ap0'-aen-Vi/^. (33) 
ar] 

The left-hand side represents the time evolution of the product {ap) for a 
photon of momentum p traveling over a geodesic. In a perfectly homoge- 
neous universe, the photon would only experience the average cosmological 
redshifting, p oc a~^, and the product (ap) would be conserved. Due to 
presence of perturbations in the universe, photons experience gravitational 
interactions and the product {ap) varies. The first term —apcj)' accounts for 
dilation, i.e. for the fact that locally, the expansion of the universe is a bit 
more advanced or delayed than the average (remember that a(l + 0) can 
be seen as the "local scale factor" ) . This means that the redshifting of the 
photon is also a bit more advanced or delayed locally. The second term ac- 
counts for the gravitational blueshifting of photons falling in gravitational 
potential wells (or redshifting of those leaving potential wells) . The energy 
e = ^^p'^ + m? can be simply replaced by p for massless photons. In that 
case we can rewrite the geodesic equation as 

)-J^ = -cj)' -fi-V^lj . (34) 

dr] 

The fact that p disappeared from the right-hand side is crucial: it shows 
that photons of different momenta traveling in the perturbed universe along 
a given geodesic all experience the same relative momentum variation. The 
consequence is that there can be no distortions of the Bose-Einstein shape 
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of /. However, photons traveling along different geodesies and in different 
directions experience different redshifting. Hence the distribution function 
acquires a dependence on one extra argument, the direction fi of propaga- 
tion (n = p/p)'. 

f{r]^x,p)= . (35) 

e^(^'-'-) - 1 

We can perform the same decomposition in terms of background and per- 
turbations as in Eq. (32). The difference is that © = ^ is now a function 
of (77, X, h). The Boltzmann equation can now be used to derive an equation 
of motion for 0(77, x,n). 

If we work in Fourier space, we can derive the equation of motion for the 
function 6(7^, ^,n). Because of the statistical isotropy of the FL universe, 
this equation does not depend explicitly on k nor n, since there is no pre- 
ferred direction: it depends only on the direction of propagation relatively 
to the considered wavenumber , i.e. on the product {k -h). Hence the equa- 
tion of motion can be written in terms of k and of the angle such that 
(k-n) = k cos 0. The initial conditions for O do depend on the wavevector k 
(since each mode gets random initial conditions) , but they also depend only 
on rather than n, for a reason that will become clear in the paragraphs 
below. Hence we can entirely eliminate h from the problem, and solve the 
equation of motion for Q{r],k,0). 

Finally, we can expand the temperature anisotropy with respect to 
using a Legendre transformation: 

e(7y,^,^) = ^(-i)H2/ + 1)9^77, ^)PKcos^) . (36) 

Here the 6^'s are the temperature anisotropy multipoles, and Pi the Leg- 
endre polynomials. It can be shown that the monopole 60 is related to the 
photon density fluctuation in a given point, the dipole 61 to its veloc- 
ity divergence and the quadrupole 62 to its anisotropic stress a^. The 
Boltzmann equation can be written as an infinite hierarchy of equations of 
motion for the coupled multipoles 

Actually, the equation of motion for 6 takes a striking form when writ- 
ten in real space, before Fourier and Legendre expansions: 

+ n • Ve - 0' + n • Vi/^ = -r (6 - Go - n • '^Te) . (37) 

We recall that F is the conformal Thomson scattering rate, 60 the temper- 
ature monopole (i.e. the average of 0(7^, x,n) over all directions n), and 
the bulk velocity of electrons, equal to that of baryons due to tight Coulomb 
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interactions. We recall that the variable is defined as the divergence of 

Vi, = Ve. 

This equation is illuminating, since it shows that at early times, when 
photons, electrons and baryons form a tightly-coupled fluid, the fact that 
r is huge forces 6 to evolve in such way that the parenthesis on the right- 
hand side vanishes. When this is the case, 6 can only have two non-zero 
component: a monopole 60, and a dipole equal to n • '^Jb. The condition on 
the dipole can be written equivalently as 0^ = ^b- 

This should remind us of the discussion of Sec. 1.3, when we noticed 
that a perfect fluid can be described in terms of 6x^ and Ox only, with a 
vanishing anisotropic stress cr^. Here, we see concretely how this conclusion 
emerges from the Boltzmann equation in the tightly-coupled regime: O2 
(related to cr^) and all higher multipoles vanish; the photon perturbations 
can be described in terms of only two independent variables 60 and 61 (or 

and Oj), like in a perfect fluid. 

The fact that in the tightly-coupled limit the system is driven towards 
0^ = 0\y simply shows that the interaction imposes a common bulk velocity 
to photons, baryons and electrons, as it should be the case in any tightly- 
coupled fluid (note that we are referring to bulk velocities, not to individual 
velocities of particles, which are still equal to c for interacting photons). 

In the tightly-coupled regime, the dipole component of 6(77, x, h) is given 
by h ' v\), with ^b = V • iTb- Hence in Fourier space this component reads 
i{h • k)k~^0]j = i{cosO)k~^0]j. This justifies the fact that initial conditions 
for 6 in Fourier space depend on 0, but not on the two degrees of freedom of 
h. As mentioned above, this property is preserved by the isotropic equations 
of motion, so that at all times we can study Q as a function of the arguments 
(77, A:, 0) instead of (77, A:, h). 



2.3. Temperature anisotropy in a given direction 

The map of temperature anisotropics that we observe today {rj = rjo) in 
our location of the universe {x = o with a proper choice of origin) when 
looking in a direction h is represented mathematically by 

^(n)=e(77,o,-n) (38) 

(since in a direction n, we see photons traveling towards —ft). Our scope is 
now to relate this quantity to perturbations on the point of the last scat- 
tering surface seen in the same direction n. This can be done by integrating 
the Boltzmann equation along the corresponding line-of-sight. 
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A good starting point consists in computing the total derivative of the 
product e^{Q-\-ip) along the trajectory of photons between the last scatter- 
ing surface and the observer. The reason for choosing this product rather 
than just 6 will become clear in a few lines: the derivative of this term will 
be easy to simplify, using the Boltzmann equation. 

The total derivative of an arbitrary function of {rj^x^fi) along the 
trajectory of photons going in a direction fi reads 

d ^ ^, dxi dT drii dT 

drj ^ ^ df] dxi df] dni 

If T is of order one in perturbations, the first two terms on the right-hand 
side are also of order one. Instead the last term is of order two, since ^ 
is of order one (this is clear from the fact that in an unperturbed universe, 
photons would travel in straight line with ^ = 0). Hence we can drop this 
term in first-order perturbation theory, and for ^ we only need to keep the 
zero-th order contribution, that can be computed assuming a homogeneous 
universe: 

^ = .. (40) 

The previous relation is just telling that photons are traveling in the direc- 
tion n and at the velocity of light: hence, in units where c = 1, dx^ = drf . 
In summary, the total derivative of T is given at first order by 

J-fry, x,n)=T' + n'VT . (41) 

drj 

Let us now replace the generic function F by 

^(ry, f , n) = e-"^^) (6(7?, f , n) + ^{r], x)) . (42) 
The total derivative of this function reads 
^ [e-"(e + ^)] = e-^ (e'^^'^h- V(e + - r'e-^ie + ^) . (43) 

We now use the linearized Boltzmann equation (37) and the fact that r' = 
— r to write the result as 

± [e--(e + ^)] = -e-W\eo + i/^ + n • ^b) + e'^c/)' + . (44) 

Finally, using the definition of the visibility function g given in Sec. 2.1, we 
get 

[e--(e + ^)] = ^ (60 + + n . ^b) + e-^(0' + ^0 • (45) 
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We can integrate this relation along the line of sight, i.e. along a straight 
line seen by the observer in a given direction —n (since the photons go in the 
direction n), starting from an early time before recombination (such that 
g-r(?7ini) ^ until the present time at which photons reach the observer 
(with by definition e"^^^^^ = 1). The result reads 



(e + 7/^)|obs= / ti7y[^(eo + ^ + n-^b) + e-"(0' + 7/;O] , (46) 



where the notation |obs means "evaluated at the observer location, along 
this line of sight", i.e. at the coordinate (7/0,0, n). 

We can gain further intuition from this equation if we use the instan- 
taneous decoupling approximation, in which all photons are assumed to 
decouple precisely at the time r/dec- In this limit, we can replace the visibil- 
ity function g by the Dirac function S]j{r] — r/dec), and by the Heaviside 
function H{r] — r/dec)- Note that the approximation g{r]) = S]j{r] — r/dec) is 
correctly normalized, since the definition of g implies J dr]g{r]) = 1. In the 
instantaneous decoupling limit, Eq. (46) reads: 



where the notation |dec means "evaluated on the last scattering surface, 
along this line of sight", i.e. at the coordinate (r/dec, — ^dec^, ^) (here rdec 
is the comoving radius of the last scattering surface). Let us now give the 
interpretation of each term in this crucial equation. 

First, B|obs is the temperature anisotropy measured by the observer in 
the direction — n, while 6o|dec is the temperature anisotropy in the point of 
the last scattering surface seen in the same direction. If only these two terms 
where present, this relation would simply tell us that that the temperature 
anisotropy seen today in a given direction is equal to the intrinsic anisotropy 
in the point where the observed photons last scattered. 

Second, the term h • Vh\dec stands for the correction to this temperature 
coming from by the usual Doppler. Indeed, this correction is caused by the 
velocity of the baryon-photon fluid (that we assumed to be tightly coupled 
until r/dec) projected along the line of sight. 

Next, we expect a correction from gravitational effects. The redshift- 
ing and blueshifting of the photons traveling along gravitational potential 
fluctuations should affect the observed temperature anisotropy relatively 
to the intrinsic one. It turns out that if the gravitational potential was 
constant in time (but, of course, not in space), this effect would be conser- 
vative, and would only depend on V^lobs ~ V^ldec- This explains the second 




(B + V^)|obs = (Go + + n . ^b)) Idee + / df] {(j)' + 




(47) 
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term on the left-hand and right-hand sides. But if ijj varies in time, the 
effect is not conservative anymore: intuitively, the amount of blueshifting 
and redshifting experienced by photons traveling across a potential well do 
not compensate each other if the gravitational well gets deeper between 
the time at which the photon enters and leaves the well. This explains the 
addition term J dr]2p. A similar effect is caused by dilation effects along the 
line-of-sight, and contributes like / drjc/). 

Finally, we can drop the second term on the left-hand side, because this 
term represents only a tiny isotropic correction to the observed anisotropics. 
It is impossible to measure it with the CMB map only, because it is formally 
equivalent to a redefinition of the average temperature f, but only by a tiny 
amount of the order of 10~^T. 

Let us write once more our result, dropping this unobservable correction, 
and grouping the terms in a particular way: 

elobs = (ep + V^)|dec + n • iTbldec + / df] {(f)' ^ xP') . (48) 

^ ^ ^ ^ ^ Vdec 

sw Doppler ^ v ^ 

ISW 

The first term is conventionally called the Sachs- Wolfe (SW) term, and 
includes the intrinsic temperature term 60 and the "gravitational Doppler 
shift" term ip at one point on the last scattering surface. The second term 
is the conventional Doppler term. The last term is called the Integrated 
Sachs- Wolfe (ISW) term and contains all non- conservative gravitational 
effects occurring in a universe with non-static metric fluctuations. 

We can gain further insight on the Sachs- Wolfe term. We will try to 
find a simpler expression for (60 + '0)1 dec, that applies at least for describ- 
ing large angular patterns on CMB maps, i.e. maps smoothed over small 
scales. For instance, this expression would describe very well the map of 
the COBE satellite, which had limited angular resolution. In more precise 
terms, we wish to calculate the contribution to the Sachs- Wolfe term of 
large wavelengths, which are bigger than the Hubble radius at the time of 
recombination. 

Let us first focus on the term So | dec- We have seen that on super- 
Hubble scales, temperature anisotropics only have a monopole and a dipole 
component, related respectively to Sj and Oj. We can be more precise now. 
We know from thermodynamics that the local value of the photon density 
is proportional to the temperature to the power four. Taking the derivative 
of \ogp^ = logT^, we get = A5T/T^ where on the right-hand side the 
temperature anisotropy is averaged over all directions n: hence 5^ = 4©o. 
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Let us now focus on the term V^|dec- We have seen in Sec. 1.4 that on 
super-Hubble scales and for adiabatic initial conditions, 5^ = |^b- Also, we 
know that decoupling takes place at the beginning of matter domination, 
and we mentioned at the very end of section 1.4 that for super-Hubble scales 
and during matter domination, one has —20 = —2?/; = d^ot = = ^cdm- 
Hence, — 2-0 = ^5^. Putting all these equalities together, we conclude that 
on the last scattering surface, 

eo + = + ^ = (-^ 2 ^ + 1) = (- ^ + 1) ^ = . (49) 

Moreover, on super-Hubble scales, we can neglect the Doppler term, 
which can be shown to be important only on sub-Hubble scales. We will see 
later that the integrated Sachs- Wolfe term plays a role on large scales, but 
only a small role, because (j) and ijj are static during most of the evolution 
after decoupling. Hence we get an approximation for CMB anisotropics 
smoothed over small scales (that was first derived by Sachs and Wolfe in 
1967): 

^ lobs, large scales — ^^l^^c ~ g^7|dec • (^^) 

In this calculation, we have seen that the term -0 wins over the term ©o, 
leading to a minus sign in front of 5^ in the above relation. This means that 
an over-density on the last scattering surface {5^ > 0), corresponding to a 
potential well {tjj < 0), leads to a cold spot in the observed map (O < 0). 
Conversely, a hot spot corresponds to an under- density. Hence, due to the 
"gravitational Doppler shift" effect accounted by the term ip (often called 
the Sachs- Wolfe effect), the patterns that we observe on CMB map are 
inverted with respect to intrinsic fluctuations on the last scattering surface. 

Equation (48) (and its large-scale approximation (50)) are important 
for pedagogical purposes, but they have no practical application: indeed, 
what we wish to calculate and to compare to observations is a theoretical 
prediction for the statistical properties of CMB anisotropics. Hence, we 
need to compute at least the CMB two-point correlation function. 

2.4. Spectrum of temperature anisotropies 

Definition, The map of CMB temperature anisotropies can be expanded 
in spherical harmonics: 

ST 

-^{h) = e(r/o,o, -n) = ^a/^l/^(n) . (51) 

Im 
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Using the Legendre expansion of O introduced in Eq. (36), and some basic 
relations between Legendre polynomials and spherical harmonics, it is easy 
to express aim as a function of G/: 



aim = i-if J ^Yim{k)eiim, k) , (52) 



where we recall that hats denote unit vectors: k = k/k. In linear perturba- 
tion theory and assuming Gaussian initial conditions, both and aim are 
Gaussian random variables. Using the orthogonality relation of spherical 
harmonics and the definitions given in Sec. 1.5, we can infer the two-point 
correlation function of the aiui s as a function of the power spectrum of G/, 
or even better, of the primordial curvature power spectrum: 

~ I f dk. 



mm' 



2^2 / ^Ql{m,k)Vn{k) 



(53) 



Here 5jj, represents the Kronecker symbol. The fact that {aimCti'm') vanishes 
for I ^ V or m ^ m' comes out of the algebra, but physically, it is a 
consequence of the homogeneity of the universe — just like the fact that all 
power spectra in Fourier space are proportional to S^ik' — k). Similarily, 
the fact that the quantity between brackets is a function of / but not of m 
is a consequence of isotropy — like the fact that in Fourier space, power 
spectra are functions of k but not k. 

The quantity between brackets is usually denoted C/, and is called the 
power spectrum of temperature anisotropics in harmonic space, or the tem- 
perature harmonic power spectrum: 

Ci^^J j^ef{vo,k)Vn{k). (54) 

In a universe with linear and Gaussian perturbations, the C/'s encode all 
the information concerning the cosmological model describing our universe 
that is contained in the CMB temperature map. 

It is worth coming back on the meaning of the averaging symbols (...) 
in Eq. (53). Since the theory of cosmological perturbation is stochastic, the 
a/^'s should be seen as random numbers, and the average is meant over 
many realizations of the theory. In a sense, "a given realization" means "a 
given universe", and the average holds over many universes, all obeying to 
the same cosmological model, which is encoded in the spectrum Ci. 

However, we observe CMB anisotropics in our universe, i.e. in only one 
realization. CMB maps allow us to measure a definite value of each aim- 
Hence the squares are not expected to be equal to C^, even if we 
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postulated the right model: there should be some scattering around Ci. 
This scattering limits the possibility to find the best theory matching the 
observations. However we can reduce considerably the scattering by noticing 
that for any fixed /, the statistical distribution of |a/^p is independent of m 
(as a consequence of isotropy, and as expressed by Eq. (54)). Hence, for an 
ideal full-sky CMB experiment, the best estimator of the underlying Q's is 
the average between all observed coefficients |a/^p with fixed /, 

-l<m<l 

In a typical universe, this quantity should be closer to the underlying Ci 
than a single |a^^*p. The way to see this mathematically is to consider 
again the theoretical (stochastic) a/^'s, and to define 

^'^^ E (56) 

-l<m<l 

By performing averages in the same sense as in equation (53), one can easily 
show that for Gaussian a/^'s, 

{Ci)=Q and {(Ci-Cif) = ^Cf . (57) 

The first equality shows that the observed C/'s defined in Eq. (55) are un- 
biased estimators of the true underlying C^'s. The second equality gives 
the typical scattering between the theory and the observations. Since for 
larger / we perform an average over more values of m, the relative scatter- 
ing decreases. This is the case with true data points, which are distributed 
qualitatively like in Fig. 1. This scattering is called cosmic variance. It can 
be seen as a theoretical error: because of cosmic variance, we cannot recon- 
struct the underlying model with infinite precision, even if we have infinitely 
precise observations. Cosmic variance is large for small Ts, meaning that 
the shape of the true underlying C/'s will always be poorly known at low /. 



Line-of-sight integral in Fourier space. According to equation (54), 
for a given primordial spectrum, the shape of the CMB spectrum Ci de- 
pends on the square of the transfer function 0/ (770 , k) . We would like to 
understand this shape at least qualitatively. In real space, we did learn a 
lot on the behavior of ©(77, x, h) by using the line-of-sight integral approach 
presented in section 2.3. A similar approach can be worked out in Fourier 
and harmonic space, i.e. for the variable Qi{r],k). We do not present here 
the intermediate steps, which can be found in Ref. 2. The final result shows 
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Fig. 1. Observed values of the temperature harmonic spectrum, C^^^ , are expected to 
be scattered around the true underlying Q's. The scattering, called cosmic variance, 
decreases with increasing I. 

many similarity with its real space counterpart, Eq. (46). It can be decom- 
posed into: 

rvo 

©K^Q, k) = df] St{v^ k) ji{k{r]o - r])) , 
^T(/?,A:)=^(6o + V^) + (^rtb)Ve-^((/)' + V^0 • (58) 

sw Doppler ISW 

We see that Qi{r]o, k) is given by the convolution of spherical Bessel func- 
tions ji{x) with a function St{t]^ /c), called the temperature source function, 
which contains the usual three terms: Sachs- Wolfe, Doppler, and Integrated 
Sachs- Wolfe. Like in the previous section, we can use the instantaneous de- 
coupling approximation, integrate the Doppler term by part, and write 
Qi{r]o,k) as: 

QiiVo^k) [Qo{Vdec^k)^ilj{r]^ec^k)]ji{k{r]o-r]dec)) 
+A:"^l9b(r?dec, k) j'i{k{r]o - Vdec)) 

+ drj [cj>'{rj,k)+i;'{rj,k)]ji{k{rjo-v)) (59) 

^dec 

(note that in the second line, i.e. in the Doppler term, the prime stands 
for the derivative of the function ji{x) with respect to its argument, not 
with respect to conformal time). This approximate result can be plugged 
into Eq. (54) to obtain the final spectrum Ci. We see that each Ci can 
be decomposed into six terms: the power spectrum Cf^ of the SW term, 
coming from the first line of Eq. (59) squared, that of the Doppler term, 
coming from the second line of Eq.(59) squared, that of the Integrated 
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Sachs- Wolfe term, coming from the third hne of Eq. (59) squared, and 
finally the three cross-spectra involving each pair of terms. 

For large values of Z, the spherical Bessel functions ji{x) and j'i{x) are 
very peaked near x c^l. Hence, for the Sachs- Wolfe and Doppler contribu- 
tions to the spectrum C^, the integral over k in Eq. (54) will pick up mainly 
modes with k{riQ — r]dec) — I- This shows that the SW contribution to Ci is 
given by the product of the primordial spectrum with the squared transfer 
function {Qq -\- ip) at a given value of rj and k: 

Cfw ^ [eo(rydec, k) + ^^(r/dec, k)f Vn{k) , k ^ ^- ^ (60) 

(for simplicity, we did not write numerical factors and powers of / or /c 
in front of this expression). In other words, Cf^ depends on the power 
spectrum of the perturbation (6o + V^), evaluated at the time of decoupling, 
and for wavenumbers in the vicinity of k = 7 ^ 5 

Cf^ ^ (IGO + '0|^)(ry,fe)~(7ydec, l/iVO-Vdec)) ' (^l) 

We reached this result with mathematical arguments, but it has a very sim- 
ple geometrical interpretation, illustrated in Fig. 2. The spectrum Ci en- 
codes the correlation between structures on CMB maps seen under an angle 
= tt/L This angle subtends a given physical scale on the last scattering 
surface, namely x (ia(^dec)7 where ds,{z) is the angular diameter distance 
to objects of redshift z. Since the Sachs- Wolfe term (Oq H-V^) contributes to 
the temperature map only at the time 77 rjdec^ the spectrum Cf^ should 
depend on the power spectrum of the Sachs- Wolfe term (|6o + at that 
time, and for a wavenumber such that 

^ = ^^=^4(.dec). (62) 

The reason for which A has been divided by two is that for spherical har- 
monics, = tt/I is the angle between a maximum and a minimum, while 
for a Fourier mode the distance between a maximum and a minimum is one 
half of the wavelength, | = ^. Eq. (62) leads to 



k 



(63) 



k I 

We recall that in a fiat FL universe the angular diameter distance is given 
by 



It{z) 



/"*° df 

,{z) = a{t{z)) -, (64) 

Jt(z') « 
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Fig. 2. A multipole I refers to structures seen on the last scattering surface under an 
angle = tt/I. These structures are seeded by Fourier modes with a half wavelength 

2 ~ k • 



t{z) being the proper time at which an object seen today with a redshift 
z emitted hght. The conformal time r]{z) is defined similarly. In terms of 
conformal time, 

d^{z) = a{r]{z)) / dr] = a{r]{z)) [r]o - r]{z)] , (65) 
and for the case of a point located on the last scattering surface, 

4(^dec) = Ci{r]dec) iVO - Vdec) • (66) 

Hence, Eq. (63) can be written as 

1 _ (^0 - ^dec) .^^x 
k~ I ' ^ ' 

which is the same relation between k and / as in Eqs. (60, 61). In those 
equations, we implicitly performed a small-angle approximation. For large 
angles (small Ts), it is inaccurate to say that a given angle/mutipole cor- 
responds to a single Fourier mode on the last-scattering surface, and it is 
important to keep the spherical Bessel function of Eq. (58) and the integral 
over k of Eq.(54). In summary, Eqs. (60, 61) represent the instantaneous 
decoupling and small- angle limit of the true power spectrum Cf^. 

Using the same two limits, a similar discussion can be carried for the 
Doppler and Integrated Sachs- Wolfe power spectra. The Doppler term de- 
pends on the power spectrum of the baryon velocity divergence evaluated 
roughly at the same time and scale, 

^Doppler ^ (|6>b|^)(ry,/c)~(rydec, l/im-mec)) ^ (68) 
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while the ISW term can be written approximately in terms of the integral 



Again, for simplicity, we did not write numerical factors and powers of / 
and k in front of these expressions. 

In the next sections, we will infer the shape of the full spectrum Ci from 
that of the three power spectra appearing in Eqs. (61, 68, 69): 



Before entering into details, we can make a guess. Even if the primordial 
spectrum Vn{k) is smooth, the temperature spectrum Ci should contain 
some structure. Indeed, as illustrated by Fig. 3, any primordial over-density 



Fig. 3. Initial over-densities in the early universe propagate in the form of wavefronts. 
The distance travelled by any wavefront at a given time is given by the sound horizon 
in the photon-baryon fluid. Here, we represent initial over-densities as spherical pat- 
terns at a given scale, while in the real universe primordial over-densities result from a 
superposition of structures on all scales. 

is expected to propagate. In Fig. 3, we represented naively the primordial 
over-densities with little dots, giving rise to spherical wavefronts at later 
time. In the real universe, wavefront patterns are less visible, because pri- 
mordial over-densities result from a superposition of structures on all scales. 
However, there is always a characteristic scale in this problem: namely, the 
distance by which a wavefront travels between some time in the primordial 
universe and the time of photon decoupling. This distance, called the sound 
horizon at decoupling (isl^dec)^ obeys to 





SW 

Doppler 
ISW 



(|6o + V^|^) at (7/,fc) (r^dec, l/ivo-Vdec)) ^ 
(|6>b|^) at (7/,/c) (r^dec, l/ivo - Vdec)) , 
(|(/>^ + V^f) for ah (7^, A:) ^(7?, l/{vo-v)) • 





(70) 
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where Cg is the sound speed in the photon-baryon fluid (in units of the 
speed of hght). Two points on the last scattering surface separated by this 
distance should be partially correlated, since density waves have propa- 
gated from one point to the other. Hence, in angular space, the two-point 
correlation function of CMB anisotropics should exhibit a characteristic fea- 
ture for angular scales corresponding to the sound horizon at decoupling, 
c^s(^dec)/c^a(^dec)- Similarly, the harmonic power spectrum Ci should 
exhibit a feature at the corresponding scale, I i\ jd i\ <^a(^dec)/<^s(^dec), 
and also for all the harmonics of this scale. We will get a confirmation of 
this in the next section. 



2.5. Acoustic oscillations 

As long as electrons, baryons and photons are tightly coupled, they form an 
effective single fluid in which density waves propagate at the sound speed 

^2 _ ^Pl + ^Ph ^^^^ 

(the density and pressure of electrons is always negligible with respect to 
that of photons). The density fluctuation Sx of each species x = 7,6 can 
be inferred from the local value of the equilibrium temperature. The fact 
that pb and p^ oc implies Sj = |(5b, and tight coupling imposes 

^7 = ^b, as we already saw in Sec. 2.2. We can simplify the expression of 
the sound speed, using also the fact that \6p\^\ <^ \Spj\. The result reads 

cf = -777^ , it = - — oc a . (72) 

It is possible to derive a simple equation of motion for the photon temper- 
ature fluctuation Go(r^,x) in the tightly-coupled regime: 

e^' + + k'clQo = -y^ + ^4>' + 4>" . (73) 



This equation follows from the combination of the continuity and Euler 
equations for photons and baryons. Given that R is proportional to the 
scale factor, we could replace by {a' /a)R = aHR. The second term on 
the left-hand side is a damping term, increasing with the contribution of 
baryons to the total energy of the fluid. The third term accounts for pressure 
forces in the effective fluid. The flrst term on the right-hand side accounts 
for the gravitational force, and the last two terms for dilation effects. 
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This equation would be that of a simple harmonic oscillator if R was a 
constant (no friction term, constant sound speed) and in absence of gravi- 
tational source terms. Then, the solution would be of the form 



with two constants of integration (6ini,^). We know that for adiabatic 
initial conditions and in the Newtonian gauge, photon density /temperature 
fluctuations should be constant in the super-Hubble limit, kr] <^ 1: this 
fixes the phase to (/p = 0. In the opposite limit, this solution corresponds 
to the propagation of acoustic oscillations. Actually, the limit between the 
constant and oscillatory regime is not set by the value of kr], but by that 
of kCsT]. In fact, the condition kcsf] ^ 1 is equivalent to X:^ d^, where A is 
a physical wavelength (A = 27ra/A:), and ds is the physical sound horizon, 
given in the case of a constant sound speed by: 



(assuming 77 ^ rjini)- Hence, the phase kc^r] of the cosine stands for the ratio 
27rds/X. Modes start oscillating when their wavelength becomes smaller 
than the sound horizon, and later on, the number of oscillations is given by 
the ratio between these two scales. 



Fig. 4. Different regions in (k, rf) space, corresponding to qualitatively different behav- 
iors for photon (and baryon) perturbations. 



In reality, R grows with time (and crosses one roughly around the time 



©0 = ©ini cos{kCsr] + (p) , 



(74) 




(75) 
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of decoupling). In addition, the gravitational source terms in Eq. (73) can 
play a role in some regimes. Let us describe qualitatively the evolution of 
Go in the different regions in (A:, rj) space shown in Fig. 4. In this figure, 
the horizontal axis corresponds to wavenumbers k (large wavelengths are 
on the left), and the vertical axis to conformal time, flowing from top to 
bottom. The super-Hubble and sub-Hubble regions are separated by the 
solid diagonal line corresponding to A: = aH (equivalent to krj = 1 during 
radiation domination). From top to bottom, the horizontal lines correspond 
to the time of equality between radiation and matter, to the time of pho- 
ton decoupling, and to the time today. The upper dashed line separates 
wavelengths bigger /smaller than the sound horizon in the baryon-photon 
fluid before decoupling (after decoupling, this notion does not make sense 
anymore). As we have just seen, this limit corresponds to A = dg, or (up to 
a factor 27r) to A: = {a/d^). At early times, Cg = l/\/3, and this condition 
reads krj = \/3. Just before decoupling, R becomes large, Cg goes to zero, 
and the comoving sound horizon {d^/a) becomes asymptotically constant, 
explaining the shape of the upper dashed line. Finally, the lower dashed 
line separates wavelengths bigger /smaller than the diffusion length defined 
in section 2.1: this line corresponds to k = l/ra, where ra is given in first 
approximation by Eq. (29). 

The evolution of 6o in the super-Hubble region is trivial: as long as 
krj <^ 1 — and a fortiori kc^rj <C 1 — the fluctuation Gq is frozen, and 
remains approximately equal to its initial value. 

The region marked with a (I) in the figure corresponds to modes that are 
crossing the sound horizon before decoupling. This is precisely the region in 
which gravitational source terms are important. They shift the zero point 
of oscillations, and boost their amplitude (due to gravitational forces and 
dilution effects). This happens during a limited amount of time, because 
the metric fluctuations quickly decay inside the sound horizon during ra- 
diation domination, making the gravitational source terms negligible. An 
approximation for the zero-point of oscillations can be found by setting Bq 
and Gq to zero in equation (73), and by keeping only the first gravitational 
term: 

^equilibrium ^ _ 1 ^ = _(1 + . (75) 

3Cg 

Since the gravitational potential is non-zero on super-Hubble scales (we 
have seen that for adiabatic initial conditions —2xp = (^tot)? the equilibrium 
point is shifted away from zero on those scales. It reaches asymptotically 
zero on sub-sound-horizon scales. 
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Region (2) corresponds to wavelengths smaller than the sound horizon 
during radiation domination. In this regime, the metric fluctuations have 
decayed, so the source term in Eq. (73) can be neglected. The friction term 
can also be neglected, because during radiation domination, R <^1. Finally 
the effective mass k'^c^ is constant in time because <C 1 implies Cg = 1/3. 
Hence we are in the simple case discussed before, and the solution is pro- 
portional to cos(A:Csr7), corresponding to stationary oscillations, symmetric 
around ©o = 0. 

Region (3) refers to wavelengths smaller than the sound horizon during 
the intermediate stage between the time of equality and that of photon 
decoupling. In this region, the metric perturbations have decayed, but R 
cannot be neglected (baryons and photons contribute to the total energy 
density with the same order of magnitude). Hence the oscillator equation 
has a non- negligible friction term (increasing with time) , and a time- varying 
effective mass (decreasing with time). The solution of the equation corre- 
sponds to damped oscillations. Physically, this damping is caused by the 
increasing inertia and decreasing pressure of the baryon-photon fluid when 
the energy density of non-relativist ic baryons takes over. 

Region refers to modes with smaller wavelength than the diffusion 
scale in the photon-baryon fluid. We have defined this scale in section 2.1. 
At early times, in the tightly-coupled limit, the mean free path of parti- 
cles in the fluid is negligible, and cosmologically interesting scales are all 
well above the diffusion length. At the approach of decoupling, the diffu- 
sion length suddenly increases, and encompasses most of sub-sound- horizon 
wavelengths. In this regime, the oscillator equaion (73) does not apply any- 
more, because we cannot describe baryons and photons in terms of a perfect 
fluid. Perturbations are then strongly damped, since diffusion tends to av- 
erage out any small-scale fluctuation. 

After describing these different regions, we are ready for understanding 
the qualitative behavior of the various relevant transfer functions, eval- 
uated at the time of decoupling. Figure 5 shows the transfer functions 
0o(^dec5 k) (solid line), ip{r]^ec^ k) (upper dotted line) and 6o(77dec5 k) (mid- 
dle dotted line) as a function of log/c. For simplicity, we neglect the role 
of the anisotropic stress generated by neutrinos (and to a lesser extent by 
photons near decoupling time). Hence we can assume (p = ip dit all times. 
For scales above the sound horizon [kcsf] <^ 1), the transfer functions are 
constant: indeed we know that, on those scales and in the Newtonian gauge, 
density and metric fluctuations are frozen. Adiabatic initial conditions im- 
pose —2ip = Stot — Sh — f — 36o. The opposite sign of 6o and ip reflects 
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Fig. 5. Transfer functions at the time of decoupling. 



the fact that an over-density corresponds to a temperature excess and a 
gravitational potential well (and vice- versa). We remember that transfer 
functions are all normalized to 7l{k) = 1 (see section 1.5): this corresponds 
to a negative 5^ (and 0o) and to a positive ip^ like in the figure. 

Because of the decay of metric fluctuations inside the Hubble radius, 
the curve smoothly decreases and tends towards zero in the small wave- 
length limit. The behavior of 0o is more complicated. Modes which are just 
crossing the sound horizon near r] = r^dec are experiencing the boost caused 
by gravitational source terms in Eq. (73): this explains the first bump in 
the solid line. For smaller wavelengths, we see oscillatory patterns corre- 
sponding to acoustic oscillations. Smaller wavelengths crossed the sound 
horizon earlier, and had more time for oscillating before decoupling. The 
maxima observed at the time of decoupling correspond to modes that could 
experience 0.5, 1, 1.5, 2, periods of oscillations before that time. The 
zero point of oscillations follows — (1 + i?)?/^, represented on the figure with a 
dashed line: this zero point reaches zero well inside the sound horizon. The 
amplitude of the oscillations is maximal for the first oscillatory pattern, 
i.e. for modes that crossed the sound horizon very recently. The second 
oscillatory pattern is reduced by the fact that those modes stayed for a 
longer time inside the sound horizon during the matter dominated regime, 
and experienced more damping due to baryons. The third and higher os- 
cillatory patterns are reduced even more by diffusion damping just before 
photon decoupling. Temperature fluctuations on very small wavelengths are 
completely suppressed by photon diffusion. 

Figure 5 also shows qualitatively the behavior of the time derivative 
0o(^dec7^)7 which exhibits oscillations that are out of phase with respect 
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to those of Oo(?7dec7^)- This wih be important in a few paragraphs, when 
discussing the Doppler effect. 

Now that we understand quahtatively the behavior of the metric and 
photon transfer functions, we can go back to the decomposition of the 
CMB temperature spectrum Ci in three terms (Sachs- Wolfe, Doppler and 
integrated Sachs- Wolfe) discussed in section 2.5. 

Sachs- Wolfe contribution. We have seen that the Sachs- Wolfe con- 
tribution to Ci is approximately given by the power spectrum of the com- 
bination {Qo-\-ip) ai T] = 77dec7 with a correspondence between k and / given 
by Eq. (67). We know that this power spectrum is given by the product of 
the primordial spectrum Vn{k) (that we can choose to be scale-invariant in 
first approximation) by the square of the transfer function (Oq + V^)^. The 
qualitative behavior of the latter has no more secrets for us. We can pick 
up the solid and upper dotted lines in figure 5, add them up, and square 
the result. The result is shown in figure 6 as a function of log k. For modes 



with kc^T] <C 1, we see a flat plateau: there, our previous calculation of the 
Sachs- Wolfe effect would apply (this part of the curve is actually called the 
Sachs- Wolfe plateau). We then observe a series of peaks. Due to the shift 
of the zero-point of oscillations given by —{1-\-R)iIj for 6o, and hence by 
—Rip for the Sachs- Wolfe term (0o + V^), there is an asymmetry between 
the first few odd and even peaks, with odd peaks being enhanced. More- 
over the overall amplitude of the peaks is suppressed in the large k limit by 
diffusion damping. A bit of algebra would show us that the envelope of the 
peaks is given in first approximation by the function exp[— (/c//cd)^], where 
is the diffusion wavenumber, related to the diffusion comoving scale of 
Eq. (29) by ^d^d = 1. 




Fig. 6. Squared transfer functions at the time of decoupling. 
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Doppler contribution. Next, we know that the C^'s receive a second 
contribution from the Doppler effect, related to the power spectrum of 
which is equal to 0^ until baryon and photons decouple from each other. It 
turns out that the photon velocity divergence 0^ is itself related to the time 
derivative of the temperature fluctuation Bq. At a very qualitative level, we 
can infer the Doppler contribution from the shape of the transfer function 
0o(^dec7^)- This contribution is null for scales above the sound horizon, 
since in this regime there are no oscillations and no significant dynamics 
in the fluid. On smaller scales, the Doppler contribution has oscillatory 
patterns, that are out of phase with respect to those of the Sachs- Wolfe 
term. The Doppler contribution is represented schematically as a dotted 
line in figure 6. 

We can now sum up the Sachs- Wolfe and Doppler contributions. Also, 
we can transpose our results for power spectra as a function of k in terms 
of Cis as a function of We have already seen in section 2.5 that there is 
a mapping betwwen the two, at least in the small- angle and instantaneous 
decoupling approximation, with a correspondence between values of k and 
/ given by Eq. (67). The result is shows in figure 7. Note that the vertical 




Fig. 7. Contribution to the Q's from the SW and Doppler terms. 

axis stands for /^Q, for the following reason. If the primordial spectrum 
was scale invariant (Vn{k) = constant) and the transfer functions were fiat 
(as it is the case for large wavelengths/small Ts), the quantity /(/ + 1)C/ 
would also be flat and independent of /. This would follow from Eq. (60) 
if we had been more carefull in keeping all the factors. It is convenient to 
plot I'^Ci or /(/ + l)Ci instead of C^, in order to display a roughly constant 
curve, just modulated by acoustic oscillations and diffusion damping. 

Diffusion damping effect. In Fig. 7, we can identify all the features 
mentioned before: the flat Sachs- Wolfe plateau, the series of oscillations 
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with enhanced odd peaks, and the exponentiahy decaying envelope of the 
peaks for large /. The envelope is now given by exp[— (///d)^], where 1^ is 
the diffusion multipole, related to the diffusion angle by /d = 7r/6>d- The 
diffusion angle is related to the diffusion scale by the usual angular 
diameter distance relation, 



At this point, we are almost done with the qualitative description of the 
CMB temperature spectrum C^ We only missed the integrated Sachs- Wolfe 
contribution, and the effect of reionization. 

Integrated Sachs- Wolfe contribution. We have seen that the in- 
tegrated Sachs- Wolfe effect is given by an integral over {i/j^ + 0') between 
photon decoupling and today. In fact, {tp^ + remains vanishingly small 
in a large part of the space (/c, r]). In Figure 8, we hatched all the regions in 
(/c, T]) space where metric fluctuations are expected to vary with time. Let 
us discuss these different regions. 



Fig. 8. Regions in (k, rf) space where metric fluctuations are expected to vary with time 
(giving rise potentiahy to an ISW effect after decouphng). 

Before photon decoupling, we know that metric fluctuations decay in- 
side the sound horizon. Instead, in the Newtonian gauge, they remain frozen 
outside the Hubble radius, excepted near times at which the equation of 
state of the universe changes: namely, at the time of equality between ra- 
diation and matter. 

Let us now discuss the variation of metric perturbations after photon 
decoupling (this is the relevant epoch for the ISW effect). Deep inside the 
matter dominated regime, one can show that metric fluctuations are static. 



Odda,{r]dec) = Ad = a(^dec)^d(^dec). 



(77) 




February 20, 2013 2:16 WSPC - Proceedings Trim Size: 9in x 6in tasi'lesgourgues 



39 

even inside the Hubble radius (at least within linear perturbation theory). 
We will justify this result in section 3.2. Hence in figure 8 there are no 
hatches during matter domination on whatever scales. Note however that 
at the beginning of matter domination, it takes some time for sub-sound- 
horizon metric fluctuations to freeze around a constant value: hence the 
hatches continue below the line corresponding to the time of equality, and 
extend till the line corresponding to photon decoupling. 

Later on, during the A (or dark energy) dominated regime, the equa- 
tion of state of the universe changes again, so metric fluctuations vary on 
all scales, like at the time of equality. A simple calculation based on Ein- 
stein equations would show that metric fluctuations are damped during this 
stage. 

In summary, contribution to the integrated Sachs- Wolfe effect can only 
come from two regions: 

• just after photon decoupling, on sub-sound-horizon scales, and 

• during A domination, on all scales. 

These two distinct contributions are usually called the Early Integrated 
Sachs- Wolfe (EISW) and Late Integrated Sachs- Wolfe (LISW) effects. One 
can show that both effects decrease with wavelength, for geometrical rea- 
sons. Hence the EISW effect is maximal for scales crossing the sound horizon 
just at the time of photon decoupling, while the LISW effect is maximal for 
the largest observable scales today. In multipole space, this means that the 
EISW effect contributes mainly to the scale of the first peak, i.e. to / ^ 200, 
while the LISW contributes mainly to the smallest multipoles / = 2, 3, 4, 
etc. The two ISW contributions are drawn on figure 9. The EISW effect 
enhances the first peak, while the LISW effect tilts the Sachs- Wolfe plateau 
even if the primordial spectrum is exactly scale invariant. 

Reionization effect. The last effect that we omitted to describe is that 
of reionization. We have seen in section 2.1 that at small redshift {z c:^ 10), 
the reionization of the universe produces a small secondary bump in the 
visibility function, corresponding physically to a small probability for CMB 
photons to rescatter at late times. This rescattering will tend to smooth out 
any temperature anisotropy pattern. Hence, reionization lowers the overall 
amplitude of the C^, but only a small amount (by approximately 15%). 
Note that the suppression of power is not uniform over the whole multipole 
range: smoothing effects cannot reach the largest observable scales (cor- 
responding to the smallest values of /). Hence the effect of reionization is 
step-like shaped, and saturates for / of the order of 40 or so (as illustrated 
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Fig. 9. ISW contribution to the temperature spectrum. 

on figure 10). 




Fig. 10. Impact of reionization on the temperature spectrum. 



2.6. Parameter dependence of the temperature spectrum 

We summarized in the last section the various physical effect contributing to 
the shape of the CMB temperature spectrum Ci. We will now recapitulate 
the effect of the various cosmological parameters on the C^'s, within the 
framework of the minimal ACDM model. 

This model assumes zero spatial curvature, and a power-law primordial 
spectrum of scalar perturbations: 

Vn{k) = A,{k/hr-' , (78) 

where is an arbitrary fixed pivot scale, is the spectrum amplitude 
at this scale, and is called the scalar tilt (the exponent is chosen to be 
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rig — 1 rather than just rig for historical reasons; with such notations, a 
scale-invariant spectrum corresponds to Us = 1). 

We recall that the Hubble parameter today, Hq, can be expressed in 
terms of a dimensionless reduced Hubble parameter h: 

Ho = 100/ikm.s-^Mpc"^ (79) 

The physical energy density of a given component x today can be expressed 
in terms of a dimensionless parameter (jJx = ^xh^'- 

pS=f2x/0°rit=f^x^=/3Wx (80) 

3( H / h)"^ 

where /3 = ^ ^ is a fixed number, with the dimension of an energy per 
volume. 

The six free parameters of the minimal ACDM model can be chosen to 

be 

{As, ns, CJb, ^m, ^A, Treio}, (81) 

where Treio is the optical depth to recombination, which is non-zero because 
of reionization in the recent universe. At first order, this is the only param- 
eter that one needs to introduce for describing reionization. For a very 
accurate description, one should introduce other parameters specifying the 
full reionization history, but these extra parameters are very difficult to 
probe experimentally. Hence, they are usually not specified. 

The above parameter basis specifies the baryon density cjb, the total 
non-relativistic (baryons + CDM) matter density cjm, and the fractional 
density of cosmological constant ^a. The photon density is implicitly as- 
sumed to match the measured value of the CMB temperature (T = 2.725 K 
implies coj ^ 2.10"^). Neutrinos are assumed to be still relativistic today, 
with an abundance relative to photons given by the prediction of the stan- 
dard neutrino decoupling model*^ . Finally, we did not include the parameter 
Hq (or h) in the parameter basis (81). Given that we are assuming a fiat uni- 
verse, h can be inferred from other parameters in (81): h = v^cj^ii70^^~^^A)- 

The parameters basis (81) is just one particular choice. Many other 
bases would be valid, for instance, {^s? ^s? ^cdm, ^o, '^reio})- The choice 
of (81) is dictated by purely pedagogical considerations: we will show that 
the shape of the CMB spectrum can easily be related to these parameters. 



^In this course, for simplicity, we do not discuss the effect of neutrinos; this effect is far 
from being neghgible, but since we consider the abundance of neutrinos as fixed, and 
their mass as irrelevant, we can explain the effect of other free parameters without taking 
neutrinos into account. Neutrino effects are described in Ref. 1. 
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Fig. 11. In the minimal ACDM model, evolution of the background density of radiation, 
matter and cosmological constant as a function of the scale factor. 

We summarize in Fig. 11 the evolution of background densities in the 
minimal ACDM model. The normalization of the radiation density pr is 
fixed by the CMB temperature and by the standard neutrino decoupling 
model. The normalization of the matter density and of the cosmological 
constant density pA is given respectively by the parameters uj^ and I^a. 
The scale factor at radiation/matter equality is given by 



Since uOr is considered as fixed, the redshift of equality is controlled by uo^ 
only. The scale factor at matter/ A equality is given by 



so it is controlled by only. The scale factor at photon decoupling has 
a small (logarithmic) dependence on cjb, that we can neglect — we will 
consider that decoupling takes place at a fixed temperature, and hence a 
fixed scale factor and redshift. 

Among the various physical effects described in the previous section, 
we can identify eight independent leading effect. We summarize them in 
Table 1, and show by which parameters they are governed. Below, we give 
more details on these effects. 

(Cl) Since all peaks in multipole space correspond to the harmonics of a 
single correlation length in real space, the scale of the peak is controlled 
(in good approximation) by a single number /peak- It depends on the 



pr = 




(82) 




(83) 
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Table 1. Independent leading effects controlling the shape of the CMB temperature power 
spectrum Ci in the minimal ACDM model. 





Effect 


Relevant quantity 


Parameter 


(CI) 


Peak scale 


^peak — 


ds\dec 
^peak <^a|dec 




(C2) 


Odd/even peak amplitude ratio 


R\dec 






(C3) 


Overall peak amplitude 


Qdec 
ao 






(C4) 


Damping enveloppe 




^dec''"d|dec 

<^a|dec 


^ CJm, ^b 


(C5) 


Global amplitude 






As 


(C6) 


Global tilt 


d log Vtz 
d log k 




Us 


(C7) 


Additional plateau tilting (LISW) 


OA 
ao 




Qa 


(C8) 


Amplitude for / > 40 only 


''"reio 







ratio of the sound horizon at decouphng by the angular diameter dis- 
tance to decouphng. The first quantity depends on the evolution prior 
to decoupling, and in particular on the expansion history and sound 
speed. Hence it depends on co^ (governing the time of equality) and 
cjb (governing Cg as a function of a). The second quantity depends on 
the expansion and geometry of the universe after decoupling, i.e. on a 
and i^o, or in our parameter basis and cjm- 
(C2) When studying the SW contribution to the C/'s, we have seen that 
the asymmetry between the amplitude of odd and even peaks depends 
on the shift of the zero-point of acoustic oscillations by a term —Rip. 
The value of the ratio R at decoupling is governed in our parameter 
basis by cJb- 

(C3) A shift in the time of radiation/matter equality affects the amplitude 
of all peaks for two reasons: it controls the duration of the intermediate 
stage between equality and decoupling, during which acoustic oscilla- 
tions are damped by baryonic effects, and the EISW. Both effects go in 
the same direction. If equality takes place earlier, there is less time for 
damping (hence all peaks are increased). The metric fluctuations are 
also less stabilized at decoupling, and the EISW is larger (hence the 
first peak is increased even more). 
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(C4) Diffusion damping near the time of recombination controls the enve- 
lope of the peaks (the function exp[— (///d)^] suppresses them, starting 
essentially from the third one). It depends on the ratio of the damping 
scale at decoupling by the angular diameter distance to decoupling. The 
first quantity depends on the Thomson scattering rate prior to decou- 
pling, and in particular on uj^ (governing the value of conformal time 
at equality) and uj\^ (governing the ionization fraction as a function of 
a). The second quantity depends on the expansion and geometry of the 
universe after decoupling, i.e. on and i^o, or in our parameter basis 
r^A and cjm- The parameter dependence of effects (CI) and (C4) could 
be thought to be similar: in fact, it is not, because the sound horizon 
and the diffusion scale depend on very different combinations of uj^ 
and cJb- 

(C5) The global amplitude of the C/'s depends trivially on that of the 
primordial spectrum, fixed by A^. 

(C6) The global slope of the C/'s depends trivially on the tilt of the pri- 
mordial spectrum, fixed by n^. 

(C7) The LISW effect tilts the Sachs- Wolfe plateau (on top of the effect of 
rig). The plateau is more lifted at small Ts if A domination is longer, 
i.e. if metric fluctuations decay during a larger amount of time. Hence 
this effect is enhanced by large values of Qa- 

(C8) The Ci amplitude is suppressed if reionization takes place early, i.e. 
if Treio is large, but without affecting the largest scales (small Z's). 

We see that in the framework of the minimal ACDM model, six param- 
eters control eight distinct physical effects with different impacts on the 
C/'s. This suggests that an accurate enough measurement of the tempera- 
ture spectrum is sufficient for fixing the six parameters of the cosmological 
model describing our universe, at least if the data is consistent with ACDM. 
This conclusion is roughly correct, but must be refined with some words of 
caution. Indeed, we remember that for small Ts, cosmic variance is large, so 
that the average Ci 's of the "true model" describing our universe cannot be 
measured precisely, even in the case of an ideal experiment. However, two 
of the previous effects (Cl) - (C8) can only affect the smallest multipoles: 

• effect (C7) affects only the Sachs- Wolfe plateau, 

• a combination of effects (C5) and (C8), corresponding to the product 
g-2rreio^^^ controls the global amplitude for / 40, but a variation 
of both Treio ciud with the previous product being kept fixed would 
only affect the smallest multipoles. 
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Hence, the measurement of the CMB temperature spectrum is sufficient for 
constraining the six parameters of the ACDM model, but with a relatively 
large error bar for Qa and for a particular combination of Treio ^ind A^. In the 
next section, we will say a few words on the measurement of CMB polariza- 
tion, which allows to better constrain reionization: polarization data allow 
to remove the degeneracy between Treio and A^. Instead, the error bar on 
can only be improved by combining CMB data with other cosmological 
probes (e.g. supernovae luminosity or large scale structure data). 

2.7. A quick word on polarization 

The Bolztmann equation presented in Eq. (37) was a bit over-simplified. 
We did as if the only degree of freedom describing photons with a black- 
body spectrum was their temperature. In fact, photons are described by 
more degrees of freedom, called the Stokes parameters, involving also their 
polarization. In Eq. (37), the Thomson scattering rate was integrated over 
polarization parameters, but in reality some polarized correction terms are 
present. 

Well before decoupling, photons remain unpolarized on average. Indeed, 
we have seen in section 2.2 that, throughout the tight-coupling regime and 
in the frame comoving with the fluid (i.e. such that 6>b = 6^7 = 0), the pho- 
ton temperature is isotropic is every point. This isotropy (resulting from 
frequent interactions) implies that photons acquire no net polarization pat- 
terns when they scatter over electrons. 

Instead, at the approach of decoupling, when Thomson scattering be- 
comes inefficient, the photon temperature is no longer isotropic in the frame 
comoving with the electrons. This means that a given electron will scatter 
simultaneously some hotter photons coming from one direction, and some 
colder photons coming from another direction. What is important for polar- 
ization is the quadrupolar component of the the temperature distribution 
in each point. This component starts from zero and grows at the approach 
of decoupling. When it becomes signiflcant, photon scattering leads to a 
net linear polarization. Hence, today, CMB photons have a different polar- 
ization amplitude and orientation in each direction of the sky. 

The map of CMB temperature anisotropics is a scalar map: it can be 
represented with a one- dimensional color code. The map of CMB polar- 
ization can be represented with sticks of different size and orientation in 
different points of the map. Roughly speaking, the size and orientation of 
the sticks can be related to the magnitude and orientation of the quadrupole 
anisotropy in each point of the last scattering surface. Since these quadrupo- 
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lar patterns reflect variations of temperature in the region of the last scat- 
tering surface, it is clear that there is a non-zero correlation between tem- 
perature and polarization maps. Still, both maps contain some independent 
information. 

In general, a vector fleld can be decomposed into a gradient and a curl 
component, like the electric and magnetic flelds. Similarly, CMB polariza- 
tion maps can be expanded in two scalar maps called, by analogy, E-type 
and B-type polarization maps. 

One can show that primordial scalar perturbations can produce both 
temperature anisotropics and E-type polarization anisotropics. If the early 
universe only features Gaussian scalar perturbations on cosmological scales, 
all the information contained in CMB maps is encoded in the temperature 
power spectrum Cj^ ^ the E-type polarization power spectrum Cf^^, and 
the cross-correlation power spectrum Cf^. The later really contains ad- 
ditional independent information, because polarization patterns are only 
partially correlated with temperature patterns (in mathematical terms, 

[cr? < crcD- 

We will mention in the next section the possibility that a significant 
amount of tensor perturbations are produced in the primordial universe. 
Such perturbations can also generate B-type polarization. If this is the case, 
aU the information is encoded in Cf^, Cf ^, Cf^ and ^. We did not 
include in this list the cross-correlation spectra C^"^^ and Cf^^ , because the 
parity symmetry imposes that they should vanish (at the level of primary 
anisotropics). 

Mapping the CMB polarization is interesting because it contains addi- 
tional information with respect to temperature anisotropics. 

First, reionization has a very distinct effect on the polarization spec- 
trum. The small fraction of photons rescattered around the time of reion- 
ization acquire additional polarization that will show up in the form of a 
peak in the small-/ branch of Cf^ (see figure 12) and Cf^^ . We have seen 
in the previous section that with temperature data only, cosmic variance 
limits our ability to measure Treio, and a fortiori to constrain additional 
parameters describing the reionization history. The effect of reionization 
in the small-/ branch of Cf^^ and Cf^ is so big that despite of cosmic 
variance, Treio can be well measured using polarization data. 

So the effect of reionization, described before as (C8), is very different for 
temperature and polarization. Instead, the other effects (Cl) - (C7) have a 
qualitatively similar impact on Cf^^ and Cf^^ . Still, measuring polarization 
is interesting, because temperature anisotropics feature a particular com- 
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Fig. 12. Typical shape of the temperature and E-type polarization power spectra in a 
ACDM model without (solid lines) or with (dotted lines) reionization of the universe at 
small redshift. 

bination of the SW, Doppler and ISW effects, while polarization provides 
information on the quadrupole 02 on the last scattering surface, which is 
correlated with the SW and Doppler terms, but probing a different combi- 
nation of them. Hence the role of polarization measurements is to remove 
some of the parameter degeneracies appearing in the analysis of tempera- 
ture data. 

Note that in the full Boltzmann equation, the equation of motion of 
temperature and polarization degrees of freedom are coupled. In the last 
sections, we neglected such a coupling. In fact, the impact of polarization 
on the evolution of temperature anisotropics is very small. Hence our qual- 
itative discussion of the various effects affecting CMB temperature remains 
valid. 

2.8. A quick word on tensors 

We have seen in section 1.1 that tensor modes are related to the traceless 
divergenceless components of the metric and stress-energy tensor, obeying 
to 

STu = and Vj, ^ diSTij = . (84) 

i i 

Since these relations impose four constraints on the six degrees of freedom 
of the 3x3 symmetric matrix STij (or 5gij)^ there are two tensor degrees 
of freedom in g^j^ and T^^^. 

The two degrees of freedom in Sgij are the two degrees of polarization of 
gravitational waves, that can propagate even in the vacuum. In the presence 
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of matter with a non-diagonal stress tensor STij, gravitational waves can 
be seeded by the tensor components STij. 

If tensor perturbations are sufficiently large at recombination or later, 
they can generate CMB temperature and polarization anisotropics: one can 
show that they interact with CMB photons, and produce effects similar to 
the SW and ISW ones. As mentioned in the previous section, gravitational 
waves can even seed B-type polarization. 

The stress-energy tensor is diagonal for perfect fluids, and negligible 
for pressureless components like CDM. Hence, during radiation and matter 
domination, tensor perturbations can only be seeded by decoupled neutri- 
nos and/or decoupled photons. However, neutrino and photon tensor modes 
are far too small for generating a detectable amount of CMB anisotropics. 
However, in the early universe, other mechanisms may produce gravita- 
tional waves on scales sufficiently large to be observable in the CMB. The 
most famous one is related to inflation. During inflation, quantum fluctu- 
ations of the metric can excite primordial tensor perturbations at a level 
that will produce a detectable signature in CMB maps. The amplitude of 
this signal is directly proportional to the energy scale of inflation. If this 
scale is large enough, tensors can contribute to Cf^^ , but only on small /'s, 
because gravitational waves decay quickly inside the Hubble radius. Hence, 
one way to detect tensors would be through a small distortion of the Ci^^ 
spectrum shape for / < 100 (i.e. for scales that are equal or larger than the 
Hubble radius at the time of decoupling). 

However, if the tensor signal is small with respect to the error bars 
associated to cosmic variance, it will remain undetectable. In that case, 
further sensitivity could be obtained by measuring C^^: in absence of ten- 
sors, would vanish, so even if the tensor amplitude is very small it can 
still dominate the C^^ signal. Unfortunately, this is true only up to some 
extent, because secondary anisotropics (in particular, those generated by 
weak lensing) produce a non-zero C^^ that could mask the primary tensor 
anisotropy spectrum. 

It is very challenging for CMB experiments to reach the sensitivity level 
required for detecting a C^^ signal, even for that generated by weak lens- 
ing. Current limits on the tensor primordial spectrum (and on the energy 
scale of inflation) mainly come from the observation of Ci^^ at small / by 
WMAP. The sensitivity of Planck to tensors will also mainly come from 
temperature. Future experiments dedicated to CMB polarization will im- 
prove the sensitivity to B-type polarization and obtain more precise bounds, 
until they reach the theoretical limit set by the lensing contamination. 
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3. Matter power spectrum 
3.1. Definition 

The total energy perturbation in the universe can be expanded as 

^^Ptot = ^Pj + + Scdm + Spr, (+ (5pde + •••) (85) 

where Sp^Q refers to possible Dark Energy (DE) perturbations, and the 
three dots for extra relics. In the minimal ACDM model, only the first four 
components are present. 

Many Large-Scale Structure (LSS) observables are related to the power 
spectrum of ^ptot, at different wavenumbers and redshifts. This is the case of 
the galaxy and of the halo correlation function, of the cluster mass function, 
of CMB lensing, of the cosmic shear spectrum, etc. 

All these observations probe the power spectrum during matter or A /DE 
domination, when photons are subdominant and Spj can be neglected. If 
neutrinos are still ultra-relativistic today, Spi, can also be neglected (in 
this course, we do not have time to discuss the impact of small neutrino 
masses, for which we refer the reader to Ref. 1). If the acceleration of the 
universe is caused by a cosmological constant, there is no term dpde- More 
generally, most Dark Energy models would predict a negligible amount of 
DE perturbations. In summary, in a wide category of cosmological scenarii 
including ACDM, we can use (5ptot — = ^Pb+^cdm- Hence, in the context 
of LSS observations, it is customary to refer to the power spectrum of the 
non-relativistic matter fluctuation defined as 



Pm Pb + Pcdm 



(86) 



which is indistinguishable from the ratio (^ptot/Pm- Only in models with 
large dark energy perturbations or modifications of Einstein gravity, the 
two quantities might be different, and special care about the definition of 
the matter power spectrum is needed. 

The matter power spectrum P{z,k) of Sm is defined like in section 1.5: 

{S^{z, k)5l^{z, k')) = 5D(k - k') P{z, k) . (87) 

Here we used the redshift as a time variable, but we could have indifferently 
used proper or conformal time. We have seen in section 1.5 that for Gaus- 
sian initial conditions and as long as perturbations are linear, the power 
spectrum at a given time can be written as the product of the primor- 
dial spectrum by the square of the relevant transfer function, in our case 
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Sm{z, k). Sticking to the same conventions as in section 1.5 and assuming a 
power-law primordial spectrum like in section 2.6, this gives 

9 2 / \ ns-i 
P{z,k) = ^A^-^J Sliz^k). (88) 

Hence, by studying qualitatively the evolution of the transfer function 
Sm{z^ /c), we will get some insight on the cosmological information encoded 
in the matter power spectrum. 

3.2. Transfer function evolution 

CDM dominated universe. In order to simplify the presentation, let us 
first assume that we leave in a ACDM universe with a negligible amount 
of baryons: <^ ^cdm and — ^cdm- In section 2.5, we wrote the mas- 
ter equation governing the evolution of photon perturbations during the 
tightly-coupled regime. Similarly, by combining the continuity and Euler 
equation of CDM perturbations, we can write here a master equation for 
(5cdm, actually valid in all regimes: 

Cdm + ^Cdm = -fcV + 3<^" + 3^<^' . (89) 

In an expanding universe, the clustering rate depends on the expansion 
rate: expansion increases distances, weakens gravitational forces, and slows 
down clustering processes. In the above equation, this is accounted by the 
second term, often called the Hubble friction term. On the right-hand side, 
the first term represents gravitational forces, and the last two terms account 
for dilation effects. 

On super-Hubble scales, we have seen that in the Newtonian gauge, 
adiabatic ICs predict constant density fluctuations Scdm- To be precise, 
and ^cdm vary on super-Hubble scales only when the total equation of state 
of the universe changes (i.e. around the time of radiation/matter equality, 
and during A domination). Instead, they remain constant on those scale 
during the radiation and matter dominated regime. 

Inside the Hubble radius, we can neglect dilation terms, and replace the 
gravitational potential term: 

'^c'dm + ^Cdm - 2 j W«)'5cdm = 0, (90) 

where ^cdm{ci) is the fraction of the critical density coming from CDM 
at a given value of time (or of the scale factor). This equation is often 
called the Meszaros equation. It can be obtained by combining Eq. (89) 
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with the (00) component of the Einstein equation (or its Poisson hmit) and 
the Friedmann equation. The careful reader might have noticed something 
suspicious: the gravitational force term k'^ijj has been eliminated in favor 
of ^cdm? while the first Einstein equation (or its Poisson limit) involves 
the total density fluctuation. Hence, shouldn't Eq. (90) feature also 5^^ 5\y 
and 5yl Actually, it turns out that the above equation is a really good 
approximation for the CDM equation of evolution in all regimes, under our 
assumption Vt\^ <C ^cdm- It applies even during radiation domination, when 
photon fluctuations are potentially large. In this course, we do not have 
enough time for justifying this point, but we refer the reader to Weinberg^'^ 
for a detailed explanation. 

During radiation domination, the Friedmann equation gives a oc ry, and 
r^cdm is much smaller than one. Hence we can neglect the last term in the 
Meszaros equation, and find that the two solutions are Scdm = constant and 
^cdm log 77. Hence CDM fluctuations grow logarithmically. During matter 
domination, a (x. rf' and f^cdm — 1, so the solutions are (^cdm cx rf' and 
^cdm 0^ Then, CDM fluctuations grow quadratically with rj. (These 
are only asymptotic solutions, but the Meszaros equation can actually be 
solved analytically at all times). During A domination, the function a{r]) 
is more complicated, and O^cdm decreases. With a bit of work, one can 
show that ^cdm grows at a smaller rate than during matter domination (i.e. 
slowlier than 77^), and that the reduction of the growth rate does not depend 
on k. 

In summary, during radiation domination, ^cdm(^5^) is constant on 
super-Hubble scales and grows logarithmically on sub-Hubble scales. A 
more precise calculation would show that up to a numerical factor of order 
one, ^cdm is given on sub-Hubble scales by ^cdm(^,^) = log(A:77). During 
matter domination, ^cdm(^, ^) is still constant on super-Hubble scales, and 
grows like rf on sub-Hubble scales. Finally, during A domination, it grows 
more slowly. These different behaviors are reported in Fig. 13. 

This simple discussion is sufficient for understanding the shape of the 
matter power spectrum at different time. Like in a cartoon. Fig. 14 shows 
this shape at four different times: at some initial time when all relevant 
modes are super-Hubble; at radiation/matter equality; at matter/ A equal- 
ity; and today. Let us comment these plots. We first need to define the 
comoving wavenumbers corresponding to wavelengths crossing the Hubble 
radius at the time of radiation/matter equality, of matter/ A equality, and 
today: 



(91) 
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Fig. 13. Qualitative evolution of the transfer function 5cdm(^5^) (normalized as usual 
to 7Z(k) = 1 at initial time) in different regimes: during radiation domination, matter 
domination, A domination, and on super/sub-Hubble scales. 

We can now review the evolution of P{k) with respect to time, fohowing 
the same steps as in Fig. 14. 

(1) At initial time, if we assume a scale-invariant spectrum with ng = 1, 
we see from Eq. (88) that P{k) oc with an amplitude given by A^. 

(2) During radiation domination, modes grow logarithmically inside the 
Hubble radius, like \og{kr]). At equality, super-Hubble modes (/Ceq^ <C 
1) are still shaped like at initial time, while sub-Hubble modes (/Ceq^ ^ 
1) have been enhanced by a factor [^cdm(^eq7 ^)/^cdm(^ini, ^)]^ — 
[log(/c?7eq)]^. The two asymptots of P{k) are then given by k~^ for 
k <C keq and /c~^[log(/c)]^ for k > /Ceq. 

(3) At the end of matter domination, when a = aA, modes still outside the 
Hubble radius keep being shaped like at initial time. This concerns all 
modes with /c <C /ca. Modes k ^ /Ceq have been amplified during mat- 
ter domination by a factor [(5cdm(^A, ^)/(^cdm(^eq, k)]'^ (^A/^eq)"^- This 
factor does not depend on k and preserves the shape of the power spec- 
trum on those scales. Finally, intermediate modes entering the Hubble 
scale during matter domination have been amplified by (t/a/^*)^, where 
T]^ is their time of Hubble crossing, given approximately by 77* = 1/k. 
Hence they have been amplified by {krij\)^. Putting all these informa- 
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Fig. 14. Shape of the matter power spectrum P(k) (log-log scale) at four different times: 
(upper left) when initial conditions are imposed (and all wavenumbers are super-Hubble); 
(upper right) at radiation/matter equality (arrows show the logarithmic growth during 
radiation domination); (lower left) at matter/ A equality (lower set of arrows show the 
growth during matter domination); (lower right) today (lower set of arrows show the 
growth during A domination). 



tions together, we see that the spectrum has three branches, scahng 
respectively hke: 

• P{k) (X for k < /ca, 

• P{k) (X k~^k^ = k for kA < k < /Ceq, 

• P{k) (X k~^{\ogk)'^ for k > /Cgq. 

(4) During A domination, (^cdm(^, ^) grows more slowly than r^^, but it still 
grows at the same rate for all sub-Hubble modes. So the shape of the 
power spectrum today is unaltered by this stage, and given by: 

• P{k) (X k~^ for k < /cq, 

• P{k) (X k~^k^ = k for /ca < /c < /Cgq, 

• P{k) (X k~^{\ogk)'^ for k > /Cgq. 

We do not enter into details for the small range of modes obeying 
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ko < k < k^. 

All this discussion was carried under the assumption of a scale- invariant 
primordial spectrum. If rig 7^ 1, the above shape should simply be rescaled 
by and the three branches of the power spectrum are given by: 

• P{k) oc k""^-^ for k < ko, 

• P{k) oc /c^" for k\ <k < /Ceq, 

• P{k) oc A:^«-4(logA:)2 for k > k^^. 

This closes the presentation of the shape of P{z^ k) in the limit 
^cdm- What remains to be seen is the impact of a non- negligible baryon 
fraction on the power spectrum. 

Baryon corrections. We must define a new important time in the 
evolution of the universe: the baryon drag time. The photon decoupling 
time was defined through the maximum of the photon visibility function 
^(77). But there are many more photons than baryons in the universe. Hence, 
for some amount of time after photon decoupling, baryons keep tracking 
photon perturbations: in other words, the baryon population decouples later 
than the photon population (seen as a whole) . 

Until the baryon drag time, we know that = |^7- After that time, 
baryons do not experience significant Thomson scattering anymore. They 
only feel gravity, and collapse in gravitational potential wells. 

We know qualitatively the behavior of the baryon transfer function 
^b(^, k) just before baryon drag, since in the CMB section we have studied 
the behavior of the photon transfer function. We know that = is con- 
stant on super-sound- horizon scales, experiences stationary oscillations on 
sub-sound-horizon scales during radiation domination, and finally damped 
oscillations on sub-sound-horizon scales during matter domination. More- 
over, we expect that at any time — and particularly at the approach of 
decoupling — the relation d\y = breaks on very small wavelengths com- 
parable to the mean free path of baryons in the imperfect baryon-photon 
fiuid, but here we will not discuss such small scales. 

The behavior of the CDM transfer function ^cdm(^5^) before baryon 
drag is still given by the Meszaros equation: deep inside the sound horizon, 
baryons are much less clustered than CDM, so that CDM is almost self- 
gravitating. 

In the left plot of Fig. 15, we sketch the qualitative behavior of the 
individual power spectrum of baryons and CDM just before baryon drag 
(Pb(^) aiid ^cdm(^) are defined as in Eq. (88), with 5^^ replaced either by 
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or ^cdm). 




Fig. 15. (Left) Shape of the baryon and CDM power spectra just before the baryon 
drag time, compared to the shape of the primordial spectrum (dashed hue). (Right) 
Shape of the matter power spectrum well after the baryon drag time, when baryons 
and CDM perturbations have reached gravitational equilibrium, compared to the shape 
of the primordial spectrum (dashed lines) and to the spectrum of a CDM-dominated 
universe (dotted line). 

After baryon drag, baryons collapse in potential gravitational wells. 
Since grows, CDM will start to feel the gravity of baryons. Finally, 
because CDM and baryons are two collisionless species feeling the same 
gravitational forces, they will equilibrate with = ^^cdm- The gravitational 
potential is then related to this common value by the (00) Einstein equation 
(or its Poisson sub-Hubble limit). This is not true however on very small 
scales, for which the baryon pressure cannot be neglected, but here we do 
not discuss such small scales. 

Hence, after a quick relaxation period, the two power spectra of baryons 
and CDM are equal to each other, Ph{k) = Pcdm(^) = P{k)' In order to 
relate this common power spectrum to the individual power spectrum of 
baryons and CDM before baryon drag, one must follow a matching process 
described in details in Einsteinstein and Hu.^ Intuitively, the matching de- 
pends on the relative weight of baryons and CDM, i.e. on the ratio l^b/^cdm- 
In the limit Q]j <C ^cdm, what we wrote at the beginning of this section ap- 
plies. In the limit l^b ^ ^cdm, the power spectrum is very suppressed with 
respect to the CDM-dominated case, with a much more negative slope on 
average, and some large oscillations corresponding to photon-baryon acous- 
tic waves before decoupling. Finally, if Q]j is a bit smaller than ftcdm but 
not negligible (as it is the case in our universe), the power spectrum de- 
parts slightly from a pure CDM one, with a smooth step-like suppression 
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plus small oscillatory patterns, which are the smoking gun in the recent 
universe of baryon-photon acoustic oscillations happening before photon 
decoupling. This is illustrated on the right plot of Fig. 15. These oscilla- 
tions are called Baryon Acoustic Oscillations (BAO). 

LSS observations can only probe the power spectrum on scales much 
smaller than the radius of the observable universe, i.e. than the current 
Hubble radius. For this reason, the first of the three branches described 
above is unobservable®. We can only measure the second and third branches, 
behaving respectively like P{k <C /Ceq) cx /c^^ and, in first approximation, 
P{k ^ /Ceq) oc [A:^^~^(log A:)^], plus the titling and superimposed oscillations 
coming from baryons. 

3.3. Parameter dependence 

The discussion presented in the previous section allows us to understand 
which effects and which parameters impact the shape of the matter power 
spectrum in the minimal (fiat) ACDM model, parametrized by 

{As, ns, CJb, ^m, ^^A, ^reio} (92) 

(see Sec. 2.6 for details on this parametrization) . We can already notice 
that the reionization optical depth is relevant for the CMB spectrum but 
not for P{k), since it only impacts the scattering rate of photons in the 
recent universe. Other effects are described in Table 2. Below, we give more 
details on these effects. 

(PI) The time of equality determines the scale k^^ of the power spectrum 
peak. More precisely, if this scale is expressed in units of /iMpc~^, 
which is the usual convention, then one can show that the scale of the 
maximum depends on both and (i-e. on our parameter basis on 
(jjm and Oa = 1 — ^m)- 

(P2) The baryon abundance (relative to CDM) is crucial at the matching 
time: a high baryon abundance leads to more suppression of the power 
spectrum for k > /cgq, and to more pronounced BAOs. 



^The first unobservable branch is actually gauge- dependent: it behaves like 
the Newtonian gauge, but not in other gauges. Note that truly observable quantities 
are always gauge-invariant. We have been a bit uncareful when saying that LSS data 
probes the power spectrum of Sm- In fact, different LSS observations probe different 
quantities, each of them being gauge-invariant. However, well inside the Hubble region, 
all these quantities coincide with each other and with the spectrum of 6m computed in 
an arbitrary gauge (up to small corrections that may become important in the future, 
but not with current experimental sensitivities). 
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Table 2. Independent leading effects controlling the shape of the matter power spec- 
trum P(k) in the minimal ACDM model. 





Effect 


Relevant quantity 


Parameter 


(PI) 


scale of the maximum 






(P2) 


slope for k ^ fceq 
and BAO amplitude 


^b/^cdm 




(P3) 


BAO scale 


r-sC^ydrag) 




(P4) 


Global amplitude 


amplitude of primordial spectrum 
and duration of AD 




(P5) 


Global tilt 


tilt of primordial spectrum 


Us 



(P3) We have seen that the scale of acoustic oscillations is set by the 
sound horizon at a given time. Since photons decouple at rydec, the 
scale of oscillations on the last scattering surface is set by ds{r]dec)i 
corresponding to the comoving scale rs{r]dec) = ds{r]dec)/ci{Vdec)- Sim- 
ilarly, since baryons decoupled at r/drag, the scale of BAOs depends 
on ds{r]drag) at baryon drag, corresponding to the comoving scale 
^s(^drag) = ds {r]drag) / ci{Vdrag) ' As explained in Sec. 2.6, the sound hori- 
zon at a given time depends on cjb and uom- In the case of ds{r]drag) the 
dependence on a;b is even stronger because the time of baryon drag 
itself depends strongly on the baryon abundance. 

(P4) The global amplitude depends of course on the primordial spectrum 
amplitude, i.e. on A^. Also, we have seen that during A domination, 
the growth rate of fluctuations is reduced with respect to matter dom- 
ination, but is independent of k for all sub-Hubble scales. Hence, 
also affects the global amplitude of the power spectrum. 

(P5) The global tilt depends of course on the primordial spectrum tilt, i.e. 
on rig. 

This discussion shows that in principle, a precise measurement of the matter 
power spectrum today (or at a given redhsift) would allow to measure 
independently cjb, ou^, ng, and Cl^ (assuming a flat ACDM universe). In 
practice, given the limited precision of current data sets, some of the effects 
described above are degenerate with each other, and matter power spectrum 
observations are mainly useful in combination with CMB observations, since 
they bring independent and complementary information. 
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Future experiments will use all the discriminating power of the matter 
power spectrum. They will perform accurate measurements of P{k, z) at 
different redshifts. This is crucial for at least two reasons: 

• First, the effect of I^a in (P5) is redshift dependent: at different red- 
shifts, fluctuations have spent more or less time during the A domi- 
nated regime, and the power spectrum amplitude has been more or less 
affected. Hence, the comparison of the amplitude at different redshifts 
allows to measure , or more generally to test the compatibility of the 
data with a cosmological constant (rather than some dynamical Dark 
Energy model). 

• Second, BAOs appear at a fixed comoving wavenumber /cbao iii Fourier 
space (related to rs{r]drag))^ but if we measure it in different LSS 
datasets at different redshifts (corresponding to different shells in real 
space) , this scale will be seen under different angles^ . By comparing the 
BAO angle at different redshifts, one can reconstruct the angular diam- 
eter distance, and therefore the expansion history at different redshifts. 
This is another way of measuring ^a or testing the A model versus DE 
models. 



4. Effects of dark matter on CMB and matter spectrum 

The existence of Dark Matter (DM) was first postulated by astronomers in 
order to account for the velocity of stars in the Milky Way (Oort, 1932), 
and that of galaxies in the Coma cluster (Zwicky, 1933). Later on, it has 
been firmly established that if general relativity applies to galactic scales 
and above, some form of DM halo must be present in most galaxies and 
clusters. 

Most of the dark matter is cold. If DM is made of collisionless and 
stable particles, the only DM property that can be probed by cosmological 
perturbations is its velocity distribution. 

At the level of homogeneous cosmology, the DM phase-space distribu- 
tion / only depends on momentum p. The number of particles in a given 
momentum bin is given by dpp^ f{p). If p^f{p) is very peaked near p = 
or the DM mass is very large, DM particles have a negligible velocity dis- 



^Here we are assuming that the BAO scale is measured in each redshift sheh transversally, 
i.e. in the direction orthogonal to the line of sight. Real experiments probe the BAO scale 
both transversally and longitudinally, so the situation is a bit more subtle than in this 
simplified discussion. 
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persion. In this case, they can be represented effectively as a pressureless 
fluid, as explained in Sec. 1.3. We recall that they are not forming a real 
fluid, due to the absence of interactions, but since their velocity dispersion 
is negligible, all particles have approximately the same bulk velocity in each 
point (described by the divergence ^cdm)- In that case the DM component 
is called Cold Dark Matter (CDM). 

If p^f{p) is not very peaked near p = or the DM mass is not very large, 
a significant fraction of DM particles may have a velocity p/m comparable 
to the speed of light, or at least to that of astrophysical objects. Typically, 
if the velocity dispersion is of the order of 10~^ (in units of c) or larger, it 
cannot be neglected when studying the CMB or LSS spectra. In that case, 
the velocity of individual particles results from two physical mechanisms: 
the dynamics imposed by gravitational forces, plus an intrinsic velocity 
distributed stochastically, according to the phase-space distribution. The 
notion of "velocity dispersion" refers to the second component, since the 
first effect tends to impose the same bulk velocity to all particles in a 
given point. DM with a large velocity dispersion cannot be described by 
fluid equations like CDM. A famous example of DM particles with a large 
velocity dispersion is neutrinos. The velocity dispersion is then fixed by the 
neutrino mass, and by their Fermi-Dirac distribution. Such DM is called 
Hot Dark Matter. Possible DM species with a smaller (but still significant) 
velocity dispersion are called Warm Dark Matter^. 

Since CDM is a very generic limit, and neutrinos are know to exist since 
the 1950's, the two most popular DM models until the 1990's were CDM and 
HDM. CDM can be motivated by many possible extensions of the standard 
model of particle physics (supersymmetry, Peccei-Quinn symmetry, etc.), 
predicting some early decoupled stable particles, with a very small velocity 
dispersion. HDM was motivated at that time by massive neutrinos with 
masses of the order of ~ 10 eV. 

These two models lead to very different predictions for the matter power 
spectrum. We have discussed before the case of CDM. A study of pertur- 
bations in a HDM (or AHDM) model would show that for k > /Ceq, the 
power spectrum P{k) would decrease with a much steeper slope. The phys- 
ical explanation is that the velocity dispersion of HDM particles prevents 
them from clustering on small scales: hence the growth rate of Sm inside 



SThe limit between HDM and WDM is usually placed in the following way: both types 
of particles have a sizable velocity dispersion, but most HDM particles become non- 
relativistic after radiation/matter equality, while most WDM particles become non- 
relativistic before that time. 



February 20, 2013 2:16 WSPC - Proceedings Trim Size: 9in x 6in tasi'lesgourgues 



60 

the Hubble radius is considerably reduced in the HDM model. 

In the 90's, the first estimates of the matter power spectrum P{k) (or of 
some related quantities) on galaxy scales and cluster scales showed clearly 
that the HDM scenario was ruled out. Today, using more accurate LSS 
data, it is firmly established that most of the DM must be cold (or eventu- 
ally, warm, but with such a small velocity dispersion that WDM would be 
indistinguishable from CDM at least on the largest observable scales). Still, 
we cannot exclude the presence of a subdominant hot or warm component. 
In fact, since we know that neutrinos have a small mass (of the order of 
~ 0.1 eV), we do expect that a few percents of DM is hot. Since the signa- 
ture of this sub-leading component has not yet been detected, the minimal 
model is stih ACDM. 

In the ACDM case, cosmological perturbations can only measure one 
DM- related parameter: its relic density, parametrized by c^cdm- This is bet- 
ter than nothing: the value of c^cdm can be used to constrain some funda- 
mental parameters governing DM production, within the framework of a 
postulated model. 

Measuring cjcdm with the CMB, Going back to Sec. 2.6, we see that 
(Jcdm can be measured "indirectly" using CMB observations: effect (C2) 
measures cjb, effect (C3) measures cc^m, and cc^cdm is given by the difference 
between these two parameters. 

There is no more "direct" effect of CDM that could be probed with 
the CMB. A priori, we could think that photons are affected by gravita- 
tional forces generated by CDM. This gravitational interaction could only 
be relevant before photon decoupling, and around the time of sound horizon 
crossing for each given wavelength (we have seen in Sec. 2.4 that well inside 
the sound horizon, photons do not feel gravitational forces anymore, be- 
cause pressure forces are much larger). But around sound horizon crossing, 
photons feel mainly the gravitational forces generated by themselves, by 
neutrinos, and to a lesser extent by baryons. In fact, in this regime, CDM 
and other components are effectively decoupled.^ 

Hence the CMB probes the CDM abundance only indirectly, as a "miss- 
ing mass", exactly like in astrophysics: it measures the total matter den- 
sity through the time of equality, the baryon density through the ratio 
between odd and even peaks, and the difference must be attributed to a 
non-relativistic component not interacting with the photons: Dark Matter. 

In this discussion, we concentrated only on effects (C2) and (C3). Indeed, 
effects (Cl) and (C4) are also sensitive to cjb and co'cdm, but only in com- 
bination with other parameters. So, unlike (C2)-h(C3), they do not bring a 
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robust and nearly model-independent evidence for DM. 

Measuring cjcdm with the matter power spectrum. The overall shape 
of P{k) for k > /Ceq brings very strong support for CDM (or eventually 
WDM with a small enough velocity dispersion). Indeed, a universe with 
only baryonic matter, or baryonic plus hot dark matter, would have a rad- 
ically different power spectrum (with a very steep slope for k > /Ceq, and 
eventually huge B AOs in the purely baryonic case) . More precisely, data on 
the matter power spectrum allows to measure the ratio uJi^/iUcdm. through 
effect (P2) of Sec. 3.3. Hence the matter power spectrum probes the CDM 
abundance independently of the CMB. This probe is more direct, since 
P{k) is really sensitive to the gravitational clustering of CDM inside the 
Hubble radius. 

Signatures of DM annihilation or decay in the CMB. Depending on 
the underlying particle physics model, DM particles could either annihilate 
or decay. Of course, the annihilation or decay rate should be very small, 
since we need a stable dark matter population. Still, a small annihilation or 
decay rate could play a role in the recent thermal history of the universe, 
and in particular in the evolution of the ionized fraction of atoms, Xe{z), to 
which the CMB is sensitive. Indeed, DM annihilation/decay products may 
include photons that could directly ionize hydrogen atoms, or particles that 
could reheat the inter-galactic medium. These mechanisms tend to increase 
Xe at a given time. Their efficiency depends on the mass of DM particles, 
on their cross-section in the case of annihilation, and on their lifetime in 
the case of decay. 

Before recombination and photon decoupling, the ionization fraction is 
anyway very close to one. In the minimal cosmological model, Xe drops 
abruptly after recombination, freezes out when particles leave chemical 
equilibrium, and increases again to one at the time of reionization, due 
to ionizing photons emitted at star formation. DM annihilation or decay 
could potentially modify this picture and play two roles: 

• at high redshift (typically 600 < z < 1100), it could slow down the drop 
of Xe after recombination, and increase its freeze-out value. This would 
appear in the CMB as a shift in the time of decoupling (impacting effect 
(Cl) and (C4) of Sec. 2.6), as a broadening of the peak of the visibility 
function g{r]) (i.e. the last scattering "surface", which is in fact more 
like a shell, would be thicker, impacting (C4)), and as an increase of the 
optical depth to recombination (impacting (C8), but not on the same 
scales as reionization). 
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• at small redshift (typically z < 40), the concentration of DM in ha- 
los could enhance its annihilation/ decay rate, and annihilation/ decay 
products could reionize atoms and increase Xg. This would be equiva- 
lent to some extra reionization: it could modify the reionization time 
or history, and impact (C8). 

The recent literature shows that for realistic models, CDM annihilation 
could play these two roles: hence, interesting bounds on DM annihilation 
can be derived from the CMB. These bounds are complementary from those 
inferred from cosmic ray data, or direct detection experiments. 

DM decay could play the second role, but the bounds on the DM lifetime 
inferred in this way from the CMB are not as strong as those from cosmic 
ray data. 

These effects had to be mentioned in this section, because they are part 
of the interplay between DM physics and cosmological perturbation theory. 
We do not discuss here bounds on DM annihilation/decay from cosmic ray 
data, since they have little to do with the evolution of linear cosmological 
perturbations. 

What else could we learn on DM from CMB or LSS data? In the 

simplest cosmological scenario, cosmological perturbations only probe the 
relic abundance of DM particles. If we are lucky, they may provide some 
hint of DM annihilation. There are other properties of DM that can be 
tested with cosmological perturbations, even if at the moment, everything 
is consistent with the minimal ACDM picture: 

• the small mass of neutrinos implies a small fraction of HDM, with 
distinct effects on the CMB and LSS power spectrum^ (effects on P{k) 
are larger: a small neutrino mass implies a step-like suppression of the 
power spectrum on small scales, by at least 4% if we assume the smallest 
possible neutrino masses that are consistent with neutrino oscillations 
experiments) . 

• a non-negligible velocity dispersion of the dominant DM component 
(like in the WDM model) or of a fraction of it (in the case of several 
DM species) would also induce a suppression of the small-scale matter 
power spectrum, at a scale and with a shape depending on the mass and 
momentum distribution of the warm species. However, such a feature 
is clearly not observed in the matter power spectrum on linear scales. 
Such a signature could only be found on smaller scales, affected by 
non-linear gravitational clustering. 

• instead of being treated as collisionless, DM could experience small in- 
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teractions with itself or with other species (photons, neutrinos baryons). 
An these possibihties require a case- by-case analysis, and usually pre- 
dict observable features similar to WDM. Both the CMB and the matter 
power spectrum can be used to put very strong bounds on a possible 
DM interactions. 

In conclusion of this section, it appears that CMB and LSS observations 
bring some crucial information concerning the abundance and the nature 
of Dark Matter. But as long as everything is compatible with the minimal 
ACDM scenario, CMB and LSS data are only telling us what DM cannot 
consist of, rather than what it is. 



5. Effects of dark energy on CMB and matter spectrum 

In this course, we assumed for simplicity that the (at least apparent) ac- 
celeration of the universe is caused by a cosmological constant. The possi- 
bility that acceleration would be caused by a dynamical component, called 
generically Dark Energy (DE), is one of the most investigated questions 
in cosmology. Of course, the acceleration could also be caused by some al- 
ternative theory of gravity, different from General Relativity on very large 
scales. Even in that case, everything could be reformulated in terms of a 
dark energy component (by moving new terms from the left to the right 
hand-side of the Einstein equation). Hence, measuring DE parameters is a 
rather generic way to look for deviations from a simple cosmological con- 
stant^ . 

Possible Dark Energy parameters. At the background level, the stress- 
energy tensor of any species compatible with the assumption of a homo- 
geneous and isotropic FL universe must be diagonal, involving only two 
function pde{ci) and pde{ci)^ related to each other by an equation of con- 
servation of energy. Hence these two functions are not independent of each 
other, at least as long as DE is not interacting with other species. Even if 
in general a DE model has no reasons to obey to an equation of state, one 
can always define the function w{a) = Pde(<^)/Pde(<^)• 
A cosmological constant corresponds to w{a) = —1, i.e. to a constant 
density pde = ~Pde^ parametrized by the value of = ^de = Pde/Pcrit 



^Note that alternatively, the observed apparent acceleration might be the consequence 
of a departure from the homogeneous Friemann-Lemaitre model, but we will not discuss 
such a possibility (which is very constrained and maybe already excluded by observa- 
tions) in this course. 
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today. 

A simple extension of tiiis model is the wCDM scenario, in which one 
assumes a linear and constant equation of state p^e = '^Pde, leading to pde ^ 
^-3(i+w) ^y. yirtue of energy conservation. This model has two independent 
DE parameters, l^de and w. 

However, w could in principle vary with time. In a purely phenomeno- 
logical approach, after confronting the wCDM model to observations, one 
can test a model with one more free parameter, accounting for the variation 
of w today. Then, there are three DE parameters, for instance ^Ide, '^(^^o) 
and ^(ao). 

At the level of perturbations, different DE models with a given back- 
ground evolution could differ through their pressure perturbation, or their 
anisotropic stress. First, let us stress that DE cannot be a perfect fluid. If 
this was the case, DE would have an adiabatic sound speed obeying to 

2 _ P' i^py , /oQ^ 
= — r = ^ = w rz r (93) 

where in the last equality we have used energy conservation. Given observa- 
tional constraints on w and w' today, the squared adiabatic sound speed of 
dark energy must be negative. If the actual squared sound speed = 
was negative, then DE perturbations would be unstable, and would ex- 
plode exponentially, in contradiction with observations. This proves that 
for DE, 7^ c^. The stability of DE perturbations imposes > 0. On 
the other hand, the absence of superluminal propagation in conventional 
physics imposes an upper bound, < 1. 

In phenomenological studies of DE models, the default assumptions con- 
cerning DE perturbations is that (i) pressure perturbations are related to 
density perturbations through = 1, and that (ii) anisotropic stress is 
negligible. This is justified by the fact that DE could be a scalar field (of- 
ten called quintessence). If this is the case, one can easily show that the 
sound speed is indeed equal to one\ while anisotropic stress only appears 
at second order in perturbations. 

A DE model with = 1 cannot lead to sizable DE perturbations. 
Indeed, = 1 implies that pressure forces in the fluid resist to gravitational 



^Actually, the quantity (5pde/<^Pde is gauge dependent. Everything becomes clearer if one 
decides to call "sound speed" the ratio ^T^de/^Pde calculated in the frame comoving with 
the DE fluid, i.e. in a gauge such that the velocity divergence ^de vanishes. If the sound 
speed of a given component derives from some fundamental principle, it should take a 
simple form precisely in this frame. Actually, it is only when this condition is imposed 
that one finds (5pde = ^Pde foi" ^ scalar field. 
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collapse. Astrophysicists like to formulate this problem in terms of the Jeans 
instability. In this language, = 1 means that the Jeans length in the fluid 
is as large as the Hubble radius. Then DE perturbations remain as small 
as in the primordial universe (of the level of 10~^), and can never leave an 
observable signature, even during DE domination. 

Only if is assumed to be very small, DE perturbations can grow and 
impact LSS observables. Indeed, if goes to zero, the Jeans length can be 
arbitrarily small, and perturbations on cosmological scales can be ampli- 
fied by the Jeans instability. In that case, one needs to postulate whether 
DE perturbations feature anisotropic stress or not. In a phenomenological 
approach, one may try to parameterize the anisotropic stress using some 
effective parameter (like a viscosity parameter). One should stress that a 
very small sound speed can be assumed, but these models sound artificial: 

<C 1 is not a generic prediction of "realistic" DE models (if any such 
model exist). 

Probing DE models with CMB and LSS data. If the DE sound speed 
is not very small, DE models can only be probed through their impact on 
the background evolution, which affects indirectly the evolution of metric 
and matter perturbations. 

Since the background expansion at late time impacts the angular diam- 
eter distance, the CMB is sensitive to DE through effects (Cl) and (C4) of 
Sec 2.6, both depending on (ia(^dec)- However, (Cl) and (C4) are also sensi- 
tive to many other quantities, so this probe is only useful in combination 
with other ones: in the best case it may remove parameter degeneracies. 

The LISW effect (C7) probes the duration of the DE dominated era, and 
could bring one bit of information on DE models. In practise, it is however 
impossible to infer precise constraints from (C7) due to the large cosmic 
variance affecting the C/'s on the scale of the Sachs- Wolfe plateau. 

We have seen in Sec 3.3 that LSS data can be useful especially if the 
matter power spectrum is measured accurately at several distinct redshifts. 
Then, both the evolution of the power spectrum amplitude (in other words, 
of the linear growth factor) and of the BAO apparent scale encode infor- 
mation on the expansion rate at several redshifts, and therefore on w or 
even w{a). This is why measurements of P(/c, z) with future ambitious pro- 
grams, like for instance the Large Synoptic Survey Telescope or the Euclid 
satellite, are expected to provide the most accurate tests of DE models. 

In summary, the CMB is not very useful for testing DE models, but can 
bring a few independent pieces of information; while LSS data at different 
redshifts is sensitive to the evolution of the background DE density. LSS 
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data should anyway be used in combination with other probes of homoge- 
neous cosmology (like supernovae luminosity) in order to get a maximum 
amount of information on DE. 

The situation that we depicted would change radically if: 

• DE would account for a significant fraction of the total energy density 
even at early times {z ^ 1), like in "early dark energy" models: then, 
DE could affect the CMB more strongly than assumed in the previous 
discussion; 

• large DE perturbations would develop as a consequence of a very small 
sound speed ^ 1: then the total power spectrum (including DE 
fluctuations) may depart significantly from the ACDM one. 

Both possibilities are still open, but no such effects are indicated by current 
data, and there are no solid theoretical motivations for our universe to fall 
in one of these two cases. 

6. Conclusions and bibliography 

We hope that this short — but very dense — introduction to cosmological 
perturbations was clear enough, and that the reader will feel like going 
further. 

Some readers might regret that, after exposing in details the effect of 
cosmological parameters on the CMB and LSS spectra (in Sees. 2.6, 3.3, 4 
and 5), we did not show concrete bounds from current experiments. There 
are three reasons for this. First, the duration of the lectures and the length 
of these notes would have exploded. Second, such information would be 
outdated very quickly by new data releases. Third, it is actually a very 
good exercise (and an efficient way to digest this course) to look at published 
bounds on cosmological parameters (inferred from CMB data, or LSS data, 
or from both) , and to figure out by oneself whether everything is consistent 
with these notes: we encourage the reader to check that the best and worst 
constrained parameter combinations can be related to the physical effects 
described here. 

We promised a bibliography in the introduction. We already mentioned 
Ref. 1, in which one can find a similar but more detailed presentation, plus 
a discussion of neutrino effects. We have cited in the core of the text three 
seminal papers (Refs. 2,3,5) providing explicit calculations on precise topics: 
the line-of-sight integral in Fourier space, ^ the justification of the Meszaros 
equation and of the fact that CDM behaves effectively as a self-gravitating 
species during radiation domination,^ and the matching between CDM and 
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baryon perturbations around the time of baryon drag, in view of computing 
analytically baryon effects in the matter power spectrum.^ 

Specialists of the field refer regularly to a concise paper of Ma & 
Bertschinger,^ containing a remarkably clear presentation of gauge issues 
and of all equations of motion in two gauges (including the Newtonian one). 

Many recent textbooks cover most topics in cosmology, including linear 
perturbations, CMB and large scale structure. We only mention in the ref- 
erence list those of S. Dodelson, V. Mukhanov^ and S. Weinberg,^ all three 
outstanding, but there exist many excellent alternatives. For a presentation 
more focused on the CMB, we recommend the thesis of W. Hu^ and the 
book of R. Durrer.^^ 
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